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Abstract. In this article we consider the anisotropic Calderon 
problem and related inverse problems. The approach is based on 
limiting Carleman weights, introduced in [13j in the Euclidean case. 
We characterize those Riemannian manifolds which admit limiting 
Carleman weights, and give a complex geometrical optics construc- 
tion for a class of such manifolds. This is used to prove uniqueness 
results for anisotropic inverse problems, via the attenuated geo- 
desic X-ray transform. Earlier results in dimension n > 3 were 
restricted to real-analytic metrics. 
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1. Introduction and main results 

1.1. Introduction. In this paper we consider the Calderon problem 
in the anisotropic case. This inverse method, also called Electrical 
Impedance Tomography (EIT), consists in determining the conductiv- 
ity of a medium by making voltage and current measurements at the 
boundary. Applications range from geophysical prospection to medical 
imaging. 
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Anisotropic conductivities depend on direction. Muscle tissue in the 
human body is an important example of an anisotropic conductor. For 
instance cardiac muscle has a conductivity of 2.3 mho in the transverse 
direction and 6.3 mho in the longitudinal direction [3]. The conductiv- 
ity in this case is represented by a positive definite, smooth, symmetric 
matrix 7 = (7-'*^(x)) in a domain Q in Euclidean space. 

If there are no sources or sinks of current in Q, the potential u in Q, 
given a voltage potential / on dQ, solves the Dirichlet problem 

^ ^1^'^) =0 inn, 



dxj \ dxk 

u = f on dQ. 

Here and throughout this article we are using Einstein's summation 
convention: repeated indices in lower and upper position are summed. 
The boundary measurements are given by the Dirichlet-to-Neumann 
map (DN map), defined by 



A f = ^jk^iy, 

' dxj 



dn 



where u = [ui, . . . , z/„) denotes the unit outer normal to dQ and u is 
the solution of the Dirichlet problem. The inverse problem is whether 
one can determine 7 by knowing A^. 

Unfortunately, A^ doesn't determine 7 uniquely. This observation is 
due to L. Tartar (see [15] for an account). Let ip : Q ^ Q he a C°° 
diffeomorphism with ip\dQ = Id where Id is the identity map. Then 



A;^ = A^ 



where 



Here ip' denotes the (matrix) differential of tp, ^ip' its transpose, and 
the dot ■ represents multiplication of matrices. 

We have then a large number of conductivities with the same DN 
map: any change of variables of Q that leaves the boundary fixed gives 
rise to a new conductivity with the same boundary measurements. The 
question is whether this is the only obstruction to unique identifiability 
of the conductivity. It is known that this is the case in two dimensions. 
The anisotropic problem can be reduced to the isotropic one by using 
isothermal coordinates (Sylvester [31]), and combining this with the 
result of Nachman [21j for isotropic conductivities gives the result for 
anisotropic conductivities with two derivatives. The regularity was 
improved by Sun and Uhlmann [30] to Lipschitz conductivities using 
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the techniques of Brown and Uhlmann [5j , and to L°° conductivities by 
Astala-Lassas-Paivarinta [2] using the work of Astala-Paivarinta [3]. 

In the case of dimension n > 3, as was pointed out in [18], this is 
a problem of geometrical nature. In this article we will focus on the 
geometric problem. 

Let (M, g) be a compact Riemannian manifold with boundary. All 
manifolds will be assumed smooth (which means and oriented. 
The Laplace-Beltrami operator associated to the metric g is given in 
local coordinates by 



dxj \ dxk 

where as usual {g-'^) is the matrix inverse of {gjk), and \g\ = det{g 
Let us consider the Dirichlet problem 

AgU = in M, u\dM = /• 

The DN map in this case is defined as the normal derivative 



'jk) 



Agf = d^u\oM = d^^-^i^k 



dM 



where z/ = i/'S^j denotes the unit outer normal to 9M, and Uk = guiJ^^ 
is the conormal. The inverse problem is to recover g from A^. 

There is a similar obstruction to uniqueness as for the conductivity. 
We have 

(LI) K^,g = Kg 

where ip is a. C°° diffeomorphism of M which is the identity on the 
boundary. As usual ip*9 denotes the pull back of the metric g by the 
diffeomorphism ip. 

In the two dimensional case there is an additional obstruction since 
the Laplace-Beltrami operator is conformally invariant. More precisely 

= l^g 

for any function c, c ^ 0. Therefore we have that for n = 2 

(1.2) Ac(^*g) = Ag 

for any smooth function c 7^ so that c\qm = 1- 

Lassas and Uhlmann [17] proved that fll.ip is the only obstruction 
to unique identifiability of the metric for real-analytic manifolds in 
dimension n > 3. In the two dimensional case they showed that (11.21) 
is the only obstruction to unique identifiability for smooth Riemannian 
surfaces. Moreover these results assume that Ag is measured only on 
an open subset of the boundary. 
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Notice that these two results don't assume any condition on the 
topology of the manifold except for connectedness. An earlier result of 
Lee-Uhlmann [18] assumed that (M, g) was strongly convex and sim- 
ply connected. The result of [T7] in dimension n > 3 was extended by 
Lassas- Taylor- Uhlmann [16] to non-compact, connected real-analytic 
manifolds with boundary. Einstein manifolds are real-analytic in the 
interior and it was conjectured in [17] that Einstein manifolds are de- 
termined, up to isometry, from the DN map. This was recently proven 
by Guillarmou-Sa Barreto [S]. 

These results on the anisotropic Calderon problem for n > 3 are 
based on the analyticity of the metric. The recovery of the metric in the 
interior of M proceeds by analytic continuation, using the knowledge 
of Taylor series of g at the boundary. Thus, these results do not give 
information from the interior of the manifold. 

On the other hand, in the isotropic case where g is a. conformal 
multiple of the Euclidean metric, many results are available even for 
nonsmooth coefficients. These results are based on special complex 
geometrical optics solutions to elliptic equations, introduced in [32]- 
These have the form 



where ( G C" is a complex vector satisfying (■( = 0, and Tq is small as 
h tends to 0. However, complex geometrical optics solutions have not 
been available in the anisotropic case, which has been a major difficulty 
in the study of that problem. 

One of the main contributions of this paper is a complex geomet- 
rical optics construction for a class of Riemannian manifolds. This 
is based on the work of Kenig-Sjostrand- Uhlmann [T3], where more 
general complex geometrical optics solutions, of the form 



were constructed in Euclidean space. Here is a limiting Carleman 
weight. 

In this paper we characterize those Riemannian manifolds which ad- 
mit limiting Carleman weights, and also characterize all such weights 
in Euclidean space. We give a construction of complex geometrical 
optics solutions on a class of Riemannian manifolds, and we use these 
solutions to prove uniqueness results in inverse problems. The inverse 
problems considered are the recovery of an electric potential and a mag- 
netic field from boundary measurements on an admissible Riemannian 
manifold, and the determination of an admissible metric within a con- 
formal class from the DN map. Let us now state the precise results of 
this article. 




u 




LIMITING CARLEMAN WEIGHTS 



5 



1.2. Statement of results. We first recall the definition of limiting 
Carleman weights. Let h > he a small parameter, and consider the 
semiclassical Laplace-Beltrami operator Pq = —h'^Ag. If (y9 is a smooth 
real-valued function on M, consider the conjugated operator 

(1.3) Po,^ = e^/'^Poe-'^/^ 

Here it is natural to work with open manifolds (i.e. manifolds without 
boundary such that no component is compact). 

Definition 1.1. A real-valued smooth function if in an open manifold 
(M, g) is said to be a limiting Carleman weight if it has non-vanishing 
differential, and if it satisfies on T*M the Poisson bracket condition 

(1.4) {p^,P</>} = when p,^ = 0, 

where is the principal symbol, in semiclassical Weyl quantization, of 
the conjugated Laplace-Beltrami operator f 1 1.3\) . 

Our first result is a characterization of those Riemannian manifolds 
which admit limiting Carleman weights. 

Theorem 1. If{M,g) is an open manifold having a limiting Carleman 
weight, then some conformal multiple of the metric g admits a parallel 
unit vector field. For simply connected manifolds, the converse is also 
true. 

Locally, a manifold admits a parallel unit vector field if and only if 
it is isometric to the product of an Euclidean interval and another Rie- 
mannian manifold. This is an instance of the de Rham decomposition 
[24j . or is easy to prove directly (see Lemma [A. 40 . Thus, if (M, g) has a 



limiting weight ip, one can choose local coordinates in such a way that 
(p{x) = Xi and 

where c is a positive conformal factor. Conversely, any metric of this 
form admits (f{x) limiting weight. 

In the case n = 2, limiting Carleman weights in (M, g) are exactly 
the harmonic functions with non- vanishing differential (see Section 2). 
The case n > 3 is more complicated. However, for the Euclidean metric 
it is possible to determine all the limiting Carleman weights. 

Theorem 2. Let Q be an open subset of R", n > 3, and let e be the 
Euclidean metric. The limiting Carleman weights in {Q, e) are locally 
of the form 

ip{x) = acpoix - xo) + b 
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where a G R \ {0} and ipo is one of the following functions: 

arg{x,uji+iuj2), 



log\x\, -r^, HTg {e'%x + i^f) , log 





with ui, UJ2 orthogonal unit vectors, 6 G [0, 27r) and ^ G R" \ {0}. 
We use the following definition for the argument function 



Remark 1.2. The possible weights are all real analytic functions. 
Some comments are made at the end of section [3] about the global 
aspect of Theorem [2l 

Let us now introduce the class of manifolds which admit limiting 
Carleman weights and for which we can prove uniqueness results in 
inverse problems. For this we need the notion of simple manifolds [[27]. 

Definition 1.3. A manifold {M,g) with boundary is simple if dM is 

strictly conwea0 , and for any point x G M the exponential map exp^ is 
a diffeomorphism from some closed neighborhood ofO in T^M onto M. 

Definition 1.4. A compact manifold with boundary {M,g), of dimen- 
sion n > 3, is admissible if it is conformal to a submanifold with 
boundary ofHx (Mo,(7o) where {Mq, g^) is a compact simple [n — 1)- 
dimensional manifold. 

Examples of admissible manifolds include the following: 

1. Bounded domains in Euclidean space, in the sphere minus a 
point, or in hyperbolic space. In the last two cases, the manifold 
is conformal to a domain in Euclidean space via stereographic 
projection. 

2. More generally, any domain in a locally conformally flat man- 
ifold is admissible, provided that the domain is appropriately 
small. Such manifolds include locally symmetric 3-dimensional 
spaces, which have parallel curvature tensor so their Cotton 
tensor vanishes (see the Appendix). 

3. Any bounded domain M in R", endowed with a metric which 
in some coordinates has the form 



Im z 



2 G C\R_. 



arg z = 2 arctan 



z\ + Re^;' 




with c > and go simple, is admissible. 



^cf. Definition EjI 
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4. The class of admissible metrics is stable under C^-small pertur- 
bations of Qq. 

The first inverse problem involves the Schrodinger operator 

Cg,g = -Ag + q, 

where g is a smooth complex valued function on {M,g). We make the 
standing assumption that is not a Dirichlet eigenvalue of Cg^q in M. 
Then the Dirichlet problem 

J Cg^qU = in M, 

I u = f ondM 

has a unique solution for any / G H^^'^^dM), and we may define the 
DN map 

^g,q ■ f ^ d^u\aM- 

Given a fixed admissible metric, one can determine the potential q from 
boundary measurements. 

Theorem 3. Let {M,g) be admissible, and let q\ and q2 be two smooth 
functions on M . If Kg^g^ = Ag^q^, then qi = q2- 

This result was known previously in dimensions n > 3 for the Eu- 
clidean metric [52] and for the hyperbolic metric [TT] . We remark that 
in the two dimensional case global uniqueness is not known even for the 
Euclidean metric. It is known for potentials coming from conductivities 
[2T] or for a generic class of potentials [2^ . 

We obtain similar uniqueness results for the Schrodinger operator in 
the presence of a magnetic field. Let A be a smooth complex valued 
1-form on M (the magnetic potential), and denote 

^g,A,q = d/'dA + q, 

where = d + iAA : C°°{M) Q^{M) and dA* is the formal adjoint 
of dA (for the sesquilinear inner product induced by the Hodge dual on 
the exterior form algebra). This reads in local coordinates 

Cg,A,gU = -\g\-'/'{d,^ +iA,){\g\'/'g^''{d,,+iAk)u) 

if A = Aj dxL 

As before, we assume throughout that is not a Dirichlet eigenvalue 
of Cg^A,q in ^1 and consider the Dirichlet problem 

r Cg^A,qU = in M, 

\ u = f ondM. 

We can define the DN map as the magnetic normal derivative 

Ag,A,q ■ f ^ dAu{u)\9M- 
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This map is invariant under gauge transformations of the magnetic 
potentiah we have 

for any smooth function %p which vanishes on the boundary. Thus, it is 
natural that one recovers the magnetic field dA and electric potential 
q from the map ^^g,A,q- 

Theorem 4. Let (M, g) he admissible, let Ai, A2 he two smooth 1-forms 
on M and let qi, q2 be two smooth functions on M. If Kg^Ai,qi = ^g,A2,q2' 
then dAi = dA2 and qi = q2- 

This result was proved in [22] for the Euclidean metric. Our proof 
is closer to [7j which considers partial boundary measurements. See 
[25] for further references on the inverse problem for the magnetic 
Schrodinger operator in the Euclidean case. 

The next result considers the anisotropic Calderon problem. Under 
the additional condition that the metrics are in the same conformal 
class, one expects uniqueness since the only diffeomorphism that leaves 
a conformal class invariant is the identity. In dimensions n > 3 this 
was known earlier for metrics conformal to the Euclidean metric [32] . 
conformal to the hyperbolic metric pjj, and analytic metrics in the 
same conformal class [T^] (based on [TS]). 

Theorem 5. Let {M,gi) and (M, (72) be two admissible Riemannian 
manifolds in the same conformal class. If ^g^ = A^j, then gi = g2- 

This article is organized as follows. In the next two sections, we study 
limiting Carleman weights and prove in particular Theorems [T] and [2J 
In Sections 4 and 5, we prove Carleman estimates and construct com- 
plex geometrical optics solutions to the Schrodinger equation. Section 
6 deals with the proofs of Theorems [3l HI and [51 The last two sections 
are devoted to two results needed in the resolution of the anisotropic 
inverse problems. The first is the injectivity of an attenuated geodesic 
X-ray transform on simple manifolds, and the second states that the 
DN map determines the Taylor expansion of the different quantities in- 
volved at the boundary. Finally, there is an appendix containing basic 
definitions and facts in Riemannian geometry which are used in this 
article. 

Acknowledgements. C.E.K. is partly supported by NSF grant DMS- 
0456583. M.S. is supported in part by the Academy of Finland. G.U. 
would like to acknowledge partial support of NSF and a Walker Family 
Endowed Professorship. We would like to express our deepest thanks to 
Johannes Sjostrand who made substantial contributions to this paper. 
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His unpublished notes on characterizing hmiting Carleman weights in 
the Euchdean case are the basis for sections 2 and 3. In particular he 
proved that the level sets of limiting Carleman weights in the Euclidean 
case are either hyperspheres or hyperplanes (see section 3). We would 
also like to thank David Jerison for helpful discussions on limiting Car- 
leman weights, and Robin Graham for useful suggestions on conformal 
geometry. 

2. Limiting Carleman weights 

We refer the reader to the appendix for a short overview on Rie- 
mannian geometry. We use (■, ■) and | ■ | to denote the Riemannian 
inner product and norm both on the tangent and the cotangent space 
and D to denote the Levi-Civita connection. Throughout this paper 
semiclassical conventions are used; we refer to [6] for an exposition 
of this theory. The principal symbol of the conjugated semiclassical 
Laplace-Beltrami operator (11. 3p is given by 

(2.1) p^ = \i\^-\d^\^ + 2i{i,d^). 

There are two reasons to use limiting Carleman weights in the con- 
struction of complex geometrical optics solutions 

u = e-^^'^+''^\a + To) 

of the Schrodinger equation (— A^, + q)u = 0. Given a function (p, 
the construction amounts to looking for solutions of the conjugated 
equation 

Po,^v + h'^qv = 

of the form v = e~'^'^(a+ro) and then applying the usual WKB method. 
This includes solving the eikonal equation 

p^{x,dilj) = 

and a transport equation on a. Note that Pq^^ is not a self-adjoint 
operator and that the symbol is complex valued. The existence of 
a solution ip to the eikonal equation implies 

Hence using limiting Carleman weights is a way to ensure that the 
former (necessary) equality is fulfilled. The other reason lies in the 
fact that one wants the conjugated operator Po,95 to be locally solvable 
in the semiclassical sense, in order to find the remainder term rg and 
go from an approximate solution to an exact solution. This means 
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the principal symbol p^p of the conjugated operator needs to satisfy 
Hormander's local solvability condition 

{p^.Pip} < when = 0. 

Since applications of the complex geometrical optics construction to 
inverse problems require to construct solutions with both exponential 
weights e'^/'^ and e~'^^'^ (in order to cancel possible exponential be- 
haviour in the product of two solutions — see section E]) and since 
= p^, it seems natural to impose the bracket condition (II .41) . 
We now proceed to the analysis of limiting Carleman weights on 
open manifolds. We first observe that the notion of limiting Carleman 
weight relates to a conformal class of Riemannian manifolds. 

Lemma 2.1. Let {M,g) be an open Riemannian manifold, and (p a 
limiting Carleman weight. If c is a smooth positive function, then ^9 is 
a limiting Carleman weight in {M,cg). In particular, if f : {M,g) 
(M, g) is a conformal transformation, then f*(p is a limiting Carleman 
weight in (M, g). 

Proof. The claim follows from the fact that the principal symbol of the 
conjugated Laplace-Beltrami operator e'^/'^AcgC^'^/'* is p^ = c~^p^ and 
that 

Both term in the right-hand side vanish when p^p = Q. □ 

Remark 2.2. This lemma gives a way to construct limiting Carleman 
weights. If we already know a limiting Carleman weight v?, then any 
function of the form o /, where / is a conformal transformation on 
{M,g), is a limiting Carleman weight. 

In particular, from the linear Carleman weight (x,^), using the in- 
version X — i> which is a conformal transformation on R" \ {0} 
(endowed with the Euclidean metric e), we obtain another limiting 
Carleman weight on M = R" \ {0}. 

Now we know that the existence of limiting Carleman weights only 
depends on a conformal class of geometries. Let a and b denote re- 
spectively the real and imaginary parts of the principal symbol of the 
conjugated operator fll.3p . so 



p^ = a + ib. 
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They are given by the following expression^ 

(2.2) a = iep-My.p = ief -|V¥.|^ 

6 = 2(d<^,0 = 2(V<^,eV 
Note that {p^^p^p} /2i = {a,b} so that the condition fll.4p also reads 

(2.3) {a, 6} = when a = b = 0. 

We start with the computation of this Poisson bracket in terms of 
the Hessian D'^ip. 

Lemma 2.3. The Poisson bracket may be expressed as 

{a, b}ix, = 4I^V(e^ + 4DV(Vy., Vy.). 

Proof. Consider ai = and 02 = Idip]"^, so that a = ai — 02, the 
Poisson bracket is given by 

{a, b} = {ai, b} + {b, 02} = Ha^b + 11^02. 

Since 02 is a function only depending on x, and b is linear in ^, the 
second bracket is easily calculated and by (lA.ip we obtain 

H^a2 = 2Lv^(|V</^P) = 4DV(V</^, V(/p). 

It remains to compute the first Poisson bracket. The Hamiltonian flow 
generated by |ai is the cogeodesic flow, 

if: — > {x{t) , C,{t)) , x{t) a geodesic, ^(t) = x^(t), 
therefore using ( ]A.2p one has 

Ha,b = 2-b{xit),m)Lo = '^9lM<t)){x{t))l^^ 

= 4— ^(x(t))|^^^ = 4DV(i(0),i(0)). 

This finishes the proof of the equality since i;(0) = □ 
The bracket condition (12. 3p now reads 

(2.4) DV(^, X) + Z^V(V^, W^) = 

for all X e T(M) such that |X|2 = |Vv?|^ (X, Vy?) = 0. 
The situation in dimension 2 is particularly simple. 

Lemma 2.4. In the case of a 2- dimensional Riemannian manifold 
{M,g), limiting Carleman weights are exactly the harmonic functions 
with non-vanishing differential. 



The musical notation is recalled in the appendix. 
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Proof. This comes from the fact that when M has dimension 2, we 
have 

D^ifiiX, X) + D^^{V^, V^) = \Vip\^ Tr Z^V = I V</?|^ Aip 

if = |V<^P and (V<^,X) = 0. □ 

We will therefore continue our investigation in dimension n > 3. The 
expression of the Poisson bracket {a, b} suggests that it is convenient 
to work with Carleman weights which are also distance functions, in 
the sense that \Vip\ = 1, since in that case D'^{p(V(p,'V(p) = 0. One 
can always reduce to this case by using the conformal metric 

(2.5) ~g = \Vv>\'g 

since the notion of limiting Carleman weights only depends on a con- 
formal class of metrics. 

Lemma 2.5. Among distance functions on an open Riemannian man- 
ifold {M,g), limiting Carleman weights can be characterized by the fol- 
lowing equivalent properties. 

(1) The Hessian of ip vanishes identically. 

(2) The gradient of ip is a Killing field. 

(3) The gradient of ip is a parallel field. 

(4) If X G M and v is in the domain of exp^ then 

fiexp^v) = (f{x) + {Vif{x),v). 

Proof. If we assume |V<^| = 1, then {Dz'Vip,'Vip) = for any vector 
field Z and the bracket condition (12.41) is equivalent to 

D^cp{X, X) = {DxVifi, X) = when (X, Vif) = 0. 

By bilinearity and by the fact that D'^(f{'Vip, Z) = {Dz'V(f,'V(f) = 
one can actually drop the orthogonality condition, and after polariza- 
tion (12. 4p is furthermore equivalent to D'^(f = 0. The equivalence of (1) 
and (2) comes from the first equality in (lA.ip . The equivalence of (1) 
and (3) follows from the second equality in flA.ip . Finally, properties 



(1) and (4) are equivalent because of the identities 
Q 

—^{exp,{tv))l^^ = D^ipiy.v) 

and of Taylor's expansion. □ 

Remark 2.6. According to (4), in Riemannian manifolds, functions 
with null Hessian are the analogue of linear Carleman weights in the 
Euclidean setting. 
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Remark 2.7. According to Lemmas 12.31 and 12.51 when ip is both a 
hmiting Carleman weight and a distance function, the bracket {a, 6} in 
(12.31) vanishes everywhere instead of just on the characteristic variety 
a = 6 = 0. 

Proof of TheoremUi Suppose that the manifold {M,g) has a hmiting 
Carleman weight ip, then is both a distance function and a limiting 
Carleman weight in {M,g) where g is the metric (12.51) conformal to g. 
According to Lemma [2. 5 j this means that VgP> is a unit parallel field. 

Conversely, if (M, g) is a simply connected Riemannian manifold 
such that (M, eg) has a unit parallel field X, then according to Lemma 
IA.3t X is both a gradient field X = Vcg'P and a Killing field. Thanks to 
Lemmas 12.51 and 12.11 this implies that ip is a. limiting Carleman weight 
in (M, eg) and in (M, g). □ 

Remark 2.8. It is now easy to justify the local coordinate expression 
in the introduction for metrics which admit limiting Carleman weights. 
If (M, g) admits a limiting weight (p, the proof of Theorem 1 shows that 
'WgP> is a unit parallel field for the conformal metric g = e~^g where 
c = iVv^l"^. Lemma [A. 41 implies that near any point of M there exist 
local coordinates such that 



One obtains, after a translation of coordinates if necessary, that ^p{x) = 
xi and 



Conversely, if g is of this form and g = e~^g, then d/dxi = Vg(xi) 
is a unit parallel field by Lemma IA.4I Thus v^(x) = xi is a limiting 
Carleman weight. 

The purpose of the next lemma is to give several properties of limiting 
Carleman weights, in particular the fact that their level sets are totally 
umbilical hypersurface^. 

Lemma 2.9. A function (p with non-vanishing differential is a limiting 
Carleman weight if and only if \V(p\~'^'Vip is a conformal Killing field. 
In particular, if (p is a limiting Carleman weight then 

(1) the Hessian of (p is determined by the knowledge of Vip and 
Dy^pVip, that is 

L>V = Ac/ + \V^\-^{d^ ® {D^^V^)^ + (Dv<^V^)^ ® dp)) 
^cf. Definition EH 
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where A = -\Vip\~'^D'^ip(y'^,Vip) , 

(2) the trace of the Hessian is given by 

/\ip = Tr^V = {n- 2)A, 

(3) the level sets of ip in (M, g) are totally umbilical submanifolds 
with normal \S/ ip\~^'W (p , with principal curvatures equal to 

/i = -|v<^|-3DV(v<^, v<^), 

(4) the eigenvalues of the Hessian are X, k and —n with 

IVv^l • 

Proof. Let g denote the metric fl2.5p conformal to g. Suppose that (p> is 
a hmiting Carleman weight in {M,g) then according to Lemma [2.51 

is a KiUing field in {M,g). This imphes that it is a conformal Killing 
field in (M, g). Conversely, if V^v? is a conformal Killing field in (M, g), 
it is both a conformal Killing field and a unit gradient in (M, g). Evalu- 
ating the equality L^-^^g = •yg at (Vgip, Vgip) gives 7 = since by flA.ll) 

^V5v^(VgV?, VgV?) = 2(Dvg^VgV?, VgV?) = L^-^\Vg(p\'^ = 0. 

This implies that Vgip is a Killing field in (M, g) hence that is a 
limiting Carleman weight in (M, g) thanks to Lemma [231 and in (M, g) 
thanks to Lemma 12.11 

Formulas flA.3l) and flA.ip give 

^|V^|"2V^^/ = 2|V<^| (^/^ <^ - rf<^ ® ^^^^ |y^|2 ^^'^^ 

If (y9 is a limiting Carleman weight, then |Vv9|~^Vv9 is a conformal 
Killing field, and the former expression equals to •yg. It is easy to 
compute 7 by evaluating the former expression at (X, X) such that 
{Vip,X) = and \X\^ = 1, and by using ([23D 

7 = 2\Vip\-^D^ip{X,X) = -2\V(p\-^D^ip(Vip,V^). 

The combination of L^Vipi-^s/ifid = 19 and those two relations gives (1). 
Besides, one also has 
2 

7 = -div (|V<^|"V(^) 

= ^(|V^r2A(/?-2|Vvpr^DV(V</?,V(/?)) 
and the two expressions of 7 yield the trace of the Hessian as in (2). 
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The level sets of are submanifolds of M, with unit normal given 
by z/ = |Vv7|~^V(y2. The tangent vectors of ^'"^{t) are orthogonal to 
Vv9 and the second fundamental form is given by 

for X and Y satisfying (X, V^^) = {Y,Vip) = 0. Thus the principal 
curvatures of f^^it) are all equal to /i. 

A consequence of the expression of the Hessian in (1) is 

(2.6) DzV^ = \Z+ I V(/^r' ((Dv^Vy^, Z)V^ + (V(/^, Z)Ds7^Vy^) . 

Since D'^Lp = \g on kei dip fl ker^D^r^VipY , A is an eigenvalue of the 
Hessian of multiplicity at least n — 2. Because of the trace, the two 
remaining eigenvalues are opposite. Considering the orthogonal vectors 

U ~r V — 

\\/ip\ \D^^Vip\ iVy^l \D^^\/ip\ 

and using (12. 6p . one has : 

DuVip = Ku, DyVip = —nv. 

One of the two vectors m, v is nonzero, therefore k and — k are the two 
remaining eigenvalues. Note that when Vy? and D^/^^Vip are dependent 
then K = |A|. In this case the eigenvalues are A (with multiplicity n — 1) 
and —A. □ 

Remark 2.10. It is an accepted fact in geometry that generic mani- 
folds in dimension n > 3 should not admit nontrivial conformal Killing 
fields, thus would not admit limiting Carleman weights. This would be 
in striking contrast with the 2D case, where every metric is locally con- 
formally flat and admits infinitely many limiting weights (the harmonic 
functions). 

In the Euclidean case, the log weight ip = log |x| plays an interesting 
role (cf. [T3], [7J and p3]). The former computation allows us to give 
a partial answer to the question: if p > is a distance function (in the 
sense that |Vp| = 1), when is cp = logp a limiting Carleman weight? 
Indeed, one has 

2 D^p dp® dp Vp A \ ^ 

D ip = — , Dy^Vip = -, and A = — , 

p p^ p-^ p^ 

thus (f = log p is a limiting Carleman weight if and only if the metric 
g has the form 

g = pD'^p + dpi® dp. 
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li if = log p is a limiting Carleman weight, then by Remark 12.81 there 
are coordinates near any point of M so that logp = Xi, \Vip\~'^ = e^^^, 
and 



Conversely, any metric of this form admits a logarithmic Carleman 
weight ip{x) = Xi = logp with p = e^^. Therefore if logp is a lim- 
iting Carleman weight, one computes the Christoffel symbols in the 
coordinates {xi,x') to be 

= 1 and Tlj = -gij{x') if i 7^ 1 or j 7^ 1- 

Thus if we denote Db = Dg^, the curvature tensor for 6 7^ 1 and 07^! 
satisfies 

= i'^i^lc ~ 9bTl^ + rf^r}^ - r'f^r^^)5fii 

A direct computation using the special form of g shows that the last 
expression vanishes. Then the sectional curvature 

^ |X|2|F|2-(X,F)2 

vanishes whenever X = di. This implies that any manifold with non- 
vanishing sectional curvature at a point p cannot admit a limiting Car- 
leman weight of the form logp near p. In particular, the sphere or 
hyperbolic space do not admit limiting Carleman weights of this form. 



3. The Euclidean case 

So far, we know at least three examples of limiting Carleman weights 
defined on open subsets of the Euclidean space: the linear weight, the 
logarithmic weight, and the hyperbolic weight 




Note that the first is defined on while the two others are on R" \ 
{0}. The purpose of this section is to determine all possible limiting 
Carleman weights for open subsets of the Euclidean space. 
The following result is well known. 
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Lemma 3.1. The only connected totally umbilical hypersurfaces in the 
Euclidean space of dimension n > 3 are parts of either hyperplanes or 
hyperspheres. 

Proof. Let S be a connected totally umbilical hypersurface in an open 
subset Q C R", let z/ denote its unit exterior normal, and let A be the 
common value of the principal curvatures. First, let us prove that A is 
constant along S. We have 

Dxiy = XX 

for all vector fields X tangent to S. Therefore we deduce 

DxDyu - DyDxv = {LxX)Y - {LyX)X + X[X, Y] . 

=D[x,Y]f 

Since the Euclidean space is flat, we have 

DxDyu - DyDxy - Dyx,Y]y = RiX, Y)u = 

therefore we deduce {LxX)Y — (LyA)X = for all vectors tangent to 
S, which means that A is constant along S. 

We now consider the vector field V = Yl^=i ^j^xj and we have 

Dx{iy - XV) = AX - XDxV = 

for all vector fields X tangent to S. This means that u — Xx is constant 
along the hypersurface. If A = 0, the normal is constant along the 
hypersurface and S is part of a hyperplane, and if A 7^ then a = 
A"^(Aa; — u) is constant along the hypersurface, and S is a part of the 
hypersphere |x — a| = 1/|A|. □ 

This provides the additional information that the level sets of a lim- 
iting Carleman weight are parts of either hyperspheres or hyperplanes. 
Let f2 be a bounded open connected subset of R", n > 3, let ip he a. 
limiting Carleman weight on Q, and consider J = if{^l) which is an 
open interval. The curvature 

of the level sets, is a smooth function of t G J. 

We will begin by proving that ip is locally of one of the forms in 
Theorem [2], thus we will suppose that f2 is a small open neighborhood 
of a point xq and J is a small interval. First suppose that /i is identically 
0, i.e. that all level sets of (f are hyperplanes. Possibly after doing a 
rotation, we may suppose that 



(3.1) 
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If Q is small enough, this implies that we can take (t, x') as coordinates 
on Q, and that the level sets of are of the form 

Xi = {uj'{t),x') + s{t) 

where as usual the prime notation stands for x' = {x2, ■ ■ ■ , x„). We will 
begin with the relation 

(3.2) ip{{uj'{t),x') + s{t),x') =t. 
Differentiation with respect to x' gives 

(3.3) uj'dx^Lp + dx'^ = on ^~^{t) 

which expresses the fact that V^p is normal to the hyperplane ip~^{t). 

Differentiating with respect to t the relation (13.21) we get the following 
equations 

(3.4) {{u',x') + s)dx,ip = l, 

(3.5) and {{Co\x') + sf dl^^ + {{Cj\x') + -s)dx,^ = f). 

It remains to compute d^^ip in order to obtain ordinary differential 
equations for uj' and s. 

Differentiating (13.31) with respect to x' gives 

upl^ip + 2u,dl^^^ip + dl^p = 0, 2<3<n. 

Remember that we have supposed /i = 0, this implies that if is har- 
monic because of (2) in Lemma [2.91 Summing up the relations above 
and using the fact that Ayj = 0, we get 

(3.6) {u"-l)dl^^ + 2{u\dl,,^) = Q. 
Differentiating (13. 3p with respect to t we also get 

[ujjdl^^ + dl^,^^) x') + s)+ djjdx,^ = 0, 2<j<n. 
Multiplying the former relations by uj and summing up yields 

(3.7) {co'^dl^^ + {co', dl^,^)) {{u', x') + s) + {u', u;')dx,^ = 0. 
The combination of (13.61) and (13. 7p gives 

(3.8) {{d;',x') + s)dlv = -^-^^dx,^. 
Finally using (13. 4p . (13. 5p and (13.81) . we have 

-^-^^{{uj',x') + s) + {{u',x') + s)=0. 

1 + Uj' 

for all t G J and all x' in a neighborhood of Xq. 
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We end up with the following system of equations 

(3.9) \ on J. 

I 1 + 

Solving the first ordinary differential equation yields 

for some constant aj. After a rotation in the x' variable, we may 
suppose that aj = for j > 3. Because of the form (13. ip of the 
gradient of at Xq, we have 

uj2(t) = tan(at + 6), = 0, 3 < j < n. 

Injecting the former solution into the equation in s leads to the follow- 
ing differential equation 

s — 2atan(at + b)s = 
which can be integrated: if a = then s{t) = ct + d and if a 7^ then 

s{t) = ctan(at + b) + d. 
Finally, this gives two possible types of limiting Carleman weights 

1. a{x — Xq, + ^ 

2. aaxg{x — Xo,uJi + iuj2) + b 
with uJi,uJ2 unit orthogonal vectors. 

Remark 3.2. The use of complex variables simplified computations 
in [7]. Note that the second weight can be written as 

a arg z + b, 

with z = {x — Xq, ui) + i{x — Xq, U2) and a, b real numbers. 

Now we assume that yU does not vanish identically on J, and con- 
sider a subinterval / C /i~^(R^) (the study on yU~^(R1) can be done 
by considering When t G /, the level sets ^^^{t) are therefore 

hyperspheres of radii r{t) = l//i(t), whose centers we denote by a{t). 
Both functions r and a are smooth on J, and we have 

(3.10) <^{a{t) + r{t)uj) =t, V(t,cj)G/xr 

for some open subset F of the unit hypersphere. The normal to the 
hypersphere v^~^(t) is 

(3.11) u=\Vip\-^Vip{a{t) + r{t)uj). 
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Differentiating tlie identity (13.101) witli respect to t, we obtain 

{Vf{a{t) + r{t)uj),a + ruj) = 1 
wliicli togetlier witli (13. lip gives 

(3.12) {Vip,uj){{a,uj) + r) =1. 

Tlien taking cu = y/\y\ in tlie former equality, multiplying it by |yp 
and differentiating with respect to y gives 

((d, uj) + r) (V^p + rLp"ijj — r{Lp"uj, uj)uj^ + (Vy^, oj) (a + fuj) = 2uj 
and thus 

(3.13) ^"u; = J^^(a + ru). 

r 

We go back to (I3.12p and differentiate it with respect to t 
{(^"uj, a + ruj) ((d, uj) + r) + (Vy^, a + fuj) = 0. 
Using (I3.13P and (13. lip , this gives 

— - I Vv^l ({a, u) +r) \a + rtuP + (a, u) + r = 
r V ^ 

=1 

which leads to the following system of equations 

{f 
a - 2-a = 
^ ^ on I. 

r - -(lap + = 

The first equation implies that the centers of the hyperspheres ip^^{t) 
are moving along a line with fixed direction k G R" 

(3.15) a{t) = ao + kj r^{s)ds 

(the indefinite integral denotes a primitive on /). Thus we have |dp = 
|A;pr^ and the second equation reads 

rr — r^|/cp — = 0. 

It is convenient to rewrite the former equation in terms of the curvature 
fi = 1/r: 

/i/i + — /i^ = 



which is equivalent to 
(3.16) det 



/i /i 
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Before proceeding to the resolution of this differential equation, let 
us first notice that it actually holds on J . Indeed, if we differentiate it 
with respect to t, then we obtain an equation on /x 



which holds on any subinterval / of /i^^(R*), but also evidently on the 
interior of yU'^(O). This implies that the equation fl3.17p holds on J by 
continuity, hence that the equation fl3.16p holds on J by integration. 
To solve the equation (13.161) there are two cases to consider. 

Case k = 0. This is the case where the hyperspheres f~^{t) are concen- 
tric. The vectors (/i, /i) and its derivative are linearly dependent, thus 
(I3.16P is equivalent to the first order differential equation /i = c/i, with 
c constant. To be more precise, either fi vanishes identically or one can 
pick to & J such that /i(to) 7^ 0, and consider the maximal interval I 
on which /i doesn't vanish. By (13.161) the derivative of on J is zero, 
therefore /i solves the equation /i = c/i on J, with c constant. Then 
the curvature is an exponential function of t which never vanishes, thus 
I = J. The function ip is easily determined to be a logarithmic weight: 



Remark 3.3. Note that the logarithmic weight can be written as 



with z = yi + i\y'\, y = x - Xq = {yi,y'). 

Case k 0. The equation (I3.17P shows that the vectors (/i, /i) and its 
second derivative are linearly dependent, this implies that /i solves a 
second order differential equation of the form /i = c/i with c a constant. 
To see this, consider the function 



which is a priori only defined outside the set of zeros of /i. Note that 
if /i(r) = then by (13.161) one has /i(r) 7^ and therefore the zeros of 
/i are isolated. Moreover, by (I3.17P one has /i(r) = 0. We can extend 
cr as a continuous function on all of J, because as t tends to r we have 



(3.17) 




ip{x) = a log \x — Xo\ + b. 



(p = aRe(log2;) + b 



a = — 



ait) 



fl(t) — /i(r) t — T 



t - T fi{t) - fi{T) 
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The Taylor expansion at higher orders actually shows that a is of class 
. Differentiating the function cr on J \ /i~^(0) gives 

o = = 

and d = on J by continuity. 

Hence the function a is constant and this implies that /i satisfies the 
equation /i — c/i = with c constant. Depending on whether c is zero, 
negative or positive, the curvature is one of the following functions 

\k\ \k\ 
(3.18) ±\k\{t + b), —sm{at + b), ±— sinh(at + 6). 

a a 

If to = ip{xo) and /i(to) 7^ 0, we may choose J so small that /x doesn't 
vanish on J. If /i(to) = 0, then either /z vanishes identically near 0, 
and that case was covered first, or is an isolated zero of ^ by (13.161) 
with k ^ 0, and we may choose J so small that fi doesn't vanish on 
J\ {to}. In conclusion, one may assume that the curvature vanishes at 
most once on J. The corresponding expressions for the centers of the 
spheres a are 

k 1 k a 

^0 ~ 77T7 r~T' <^0 



\k\^t + b' ^ |A;|2 tan(at + 6)' 
k a 
" tanh(at + b) 

If /i(to) = then these expressions depend a priori on the connected 
component I± of J\ {to} on which /i is positive or negative. The norm 
of the vector k however is determined by the equation (I3.16P which 
holds on J and its direction k/\k\ by the normal to the hyperplane 
(p~^{to). This may be seen in the following way: consider the family of 
hyperspheres 



X — a± + 



k± a ^ -2 



a 



\k\^ smh\at + b) 



I A; 1 2 ta.nh{at + b) 

when t E I±. After expansion, we obtain 

\k\Hanh{at + b) ( 
(3.19) 2(a;-a±,fc±) = i ^ TTH ~ "~ 

and letting t tend to to = —h/a^ we get the equation of the hyperplane 

<^-Hio) 

(x — a±, k±) = 0. 

This implies that either k^ = k^ or k^ = —k^. Besides, by (13. lip and 
(I3.12p . we have {a, ui) + r = |V</9|~^ on /_|_, thus V^p) — /i|Vv9| = /i^ 
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and by letting x tend to xq, this gives 

(Vv?(xo),/c+) = /i(to)|Vv3(xo)|. 

Considering /_, we obtain the same identity with fc+ replaced by k_. 
This additional information removes the uncertainty on the sign and 
we have = k_. Moreover the components of a+ and a_ along k are 
also equal. The computations are similar in the two other cases. 

To determine explicitly the function ip, let us deal with the case 
where the curvature is sin(at + b)/a: to fix the ideas, suppose that 
aJ + b C. {—9, —6 + 27r) with 6' G R and that vanishes at to = 
(the non- vanishing case is simpler). Then we have ^ \ {to} = /- U /_|_ 
and 

a a 
X = a± — -— — — k + — — uj 

tan(ay9(xj + oj fc sm(av9(xj + oj 



for all X G '^{I±)- Thus, if we take ^ = —ak/\k\'^, we have 



and since aip + b + 9 — n G (— vr, tt), and (d, cj) + f > by ( 13.1ip we 
have 

aip{x) + 6 = arg {e'%x - a± + i^f) + vr - 6* 
on </)~^(J-|-). With the former expression of (p we have 

2|^| 

\aV<^{x)\ = —— _ on {x-a±,O = 

this shows that |x — a+p = |x — in a neighborhood of xq in the 
hyperplane (p~^{to), hence that = The former expression of (f 
therefore holds on O. 

The two other cases are similar. This finally gives three possible 
types of limiting Carleman weights 

{x-xo,0 

1. a- r-r:r + t> 

\x — Xo| 

2. a SiTg {e'%x ~ xo + iO^) + b 

2{x-xo,0 , . 1 k-Xo + ^P 
6. a arctanh — — -- + b = a log z -777 + b. 



Remark 3.4. As in Remark 13.31 these functions take a simple form 
with respect to some complex variable. Take y = x — xq, and {yi,y') 
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such that y = yi^/\C, \ + y' where y' is orthogonal to ^. Denote z = 
Hi + ^ly'l; the Carleman weights take the form 

1. a Re — \- b, 

z 

^ , e^^(z + ic) 

2. almlog — ^ - + b, 

z — ic 

z ~\~ c 

3. a Re log 1-6, 

z — c 

with a, b, c real numbers. 

This proves Theorem [2l Note that this result allows to determine 
limiting Carleman weights locally on the sphere or on the hyperbolic 
space by conformal transformation. 

Remark 3.5. While the linear weight is smooth on R", the weights 
log |x| and (x, are only well defined and smooth on R" \ {0}, 

and the function log \x + — on R" \ — ^}. The weight 

arg(a:, ci;i + iuj2) is defined on R" \ {{x,uj2) = 0, {x,uji) < 0} and the 
function arg (e*^(x + i^Y) on 

R"\ ||x + cotan^er = (a;,0|o| 
L sm 6 ) 

if 6^ G (0, vr) (with the > inequality) or 6' G (tt, 27r) (with the < inequal- 
ity) and on 

R"\{(x,O = 0, \x\ I 1^1} 
if ^ = (with the < inequality) or = vr (with the > inequality). 

Let us end this section with some comments about the global aspect 
of Theorem [2] when f2 is an open connected set. Let be a limiting 
Carleman weight on (f2, e). The function y9 is real analytic on f2: indeed 
if G f2, the function is one of the limiting Carleman weights 
calculated above in the neighourhood of xq, hence real analytic. This 
ensures that if is locally equal to one of the functions 

a(x — xo,0 + ^5 alog |x — xo| + 6, 

a rT> alog-^ T7^ + o 

\xY \x — xq — 41 

then if) is equal to this function on the whole set f2. Indeed if ^ is equal 
to a linear weight near some point Xq G then is equal to this weight 
on the whole set f2 by analytic continuation. If is equal to the function 
a log I X — XqI + 6 near Xq G then f2 cannot contain the singularity 
xo of this function. Otherwise, would be equal to this function on 
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fi\{a;o} by analytic continuation and would not blow up at xq. Thus (f 
is equal to a log \x — Xq \ +b on Q. The proof is similar if ip is one of the 
functions a{x — Xo,C,)/\x — Xo\'^ + b or a\og{x — Xq+^Y / {x — Xq — ^Y + b. 

For the two argument forms, we need some additional assumptions 
on fl. Suppose that ip is equal to aarg (e*^(x — Xq + i^)^) + b near xq, 
and that the image of the set fl hj x ^-^ {x — xq + i^Y is contained 
in a simply connected set U G C*, then by analytic continuation, ip is 
equal to 

a a.Tg^{x — Xq + + b 

where arg^; is the determination of the argument on U which coincides 
with arg(e*^2;) at zq = {xq — xq + iC)^- The point is similar for the 
function a arg(a; — xq, oji + 1002). In particular, if VL is contained in one 
of the domains of existence computed in Remark 13.51 and (p is locally 
of the corresponding argument form, then this is still true globally on 



4. Carleman estimates 

Let (M, g) be a compact Riemannian manifold with boundary. By 
dV we denote the volume form on {M,g), and by dS = ujdV the 
induced volume form on dM. The norm of a function is then given 
by 

lkllL2(Af) = ( / \u\'^dV 
^ Jm 

and the corresponding scalar product by 

{u\v) = / uvdV. 
Jm 

Similarly on the boundary, the norm and scalar products are given by 
\LHdM) = (l \f?dsY {f\hU,= [ fhdS. 

^ 'dM ^ JdM 



We write for short 



|VM||i2(M) = |||Vm|||^2(a^) = I / Wu?dV 



M 



and we denote by Hl^^{M) the semiclassical Sobolev space associated 
to the norm 

||«||ffi^,(M) = {\\u\\l^{M) + \\h^Al^M)Y- 

We assume that (M, g) is embedded in a compact manifold (A^, g) with- 
out boundary, and that is a limiting Carleman weight on ([/, g) where 
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U is open in and M U. The goal of this section is to prove the 
following Carleman estimate. 

Theorem 4.1. Let {U,g) be an open Riemannian manifold and {M,g) 
a compact Riemannian submanifold with boundary such that M (£ U. 
Suppose that (p is a limiting Carleman weight on {U,g). Let X be 
a smooth vector field on M and q a smooth function on M. There 
exist two constants C > and < Hq < 1 such that for all functions 
u G C^{M°) and all < h < Hq, one has the inequality 

(4.1) \\ehu\\Hi^jM) < Ch\\eh (A + X + q)u\\L2(M)- 

To lighten the notations we will forget the subscript L'^{M) whenever 
it is not needed. 

Proof. We first observe that the result is invariant under conformal 
change of metrics since we have 

c—[Ag + X + q)u = ( Ac-ig + cX + (c— m) 

with qc = cq — ^^Xc + c^Ag(c ~). Therefore if needed, we can 
assume the limiting Carleman weight to be a distance function by re- 
placing g by the conformal metric (12.51) . 

Our next observation is that the estimate may be perturbed by zero 
order terms since this gives rise to an error of the form (9(/i)||e'^/'^M||, 
which may be absorbed in the left-hand side if h is assumed small 
enough. Therefore we can neglect the potential q, and assume that 
g = from the start. Let us first assume that we also have X = 0. 
Then the Carleman estimate (14. ip is equivalent to the following a priori 
estimate 

(4.2) \\v\\H^jM)<Cih-'\\Po,^v\\. 

One goes from one inequality to another by taking V = e'pl^u. The 
conjugated operator is given by 

Po,^ = -/^^A - I V^p + 2(V<^, W) + /lAvJ. 
Then we have in particular 

= iPo,^v\v) + |||V(/?|t;||^ - 2{{V(fi,hVv)\v) - h{A(fiv\v) 
therefore using Cauchy-Schwarz inequality one sees that 

(4.3) \\hVvf < ||Po,^^ir + C^i||^^f • 

This means that the gradient of v may be controlled and that it suffices 
to prove the a priori estimate 

ll^ll < C2h-^\\Po^^v\\ 
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to obtain (Oil . 

We decompose P^^^p into its self-adjoint and skew-adjoint parts 

Po^^ = A + iB, A = -/i^A - |V(^|^ 5 = 2(V<^,/iV) + /iA(^ 

and we have by integration by parts 

(4.4) ||^o,^^^ir = \\Avf + \\Bvf + i{[A, B]v\v). 

A direct application of the commutator method will not be enough 
to get an a priori estimate assuming the bracket condition (12.31) . one 
needs to use convexification. This classical argument consists in taking 
a modified weight / o where / is a convex function chosen so that 
the bracket in (12. 3p becomes positive. We decompose the operator 
Pojoifi = A + iB into its self-adjoint and skew-adjoint parts, and denote 
by a and b the corresponding principal symbols. We now suppose, as 
we may according to our first observation, that (f is both a limiting 
Carleman weight and a distance function. We have 

V(/ovp) = (/'ov9)V^ 

=0 

therefore using Lemma [2.31 

{a, 6}(x, = 4(r o if o ^)2| v^r + 4(r o ^) (v^, 

^ V ' 

At the operator level, this gives 

i[A, B] = Ah{f" o if) (/' o if)^ + hB(3B + h^R 

where i? is a first order semiclassical differential operator. For the 
function /, we choose the following convex polynomial 

fis) = s + ^/, r is) = 1 + ^8, ns) = ^. 

We choose h/e < Eq < 1 with Eq small enough so that /' > | on ip{M) 
and denote ipi; = f o ip. Note that the coefficients of R, as well as /3, 
are uniformly bounded with respect to h and e. 
We finally obtain 

i{[A,B]v\v) > —\\vf - CshWBvf - C3h^\\v\\HiJv\\ 
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and using (14 .Sp 

z{[A,B]v\v) > ^(1 - C,e)\\vr - CsHEvf - Cs\\Po,^M 
Going back to (14. 4p . this gives 

(4.5) (l + C3)||Po,^.^f > \\Avf + il-Csh)\\Bvf 



+ -{1-C^e)\\v 



|2 



If we don't assume X = 0, the conjugated operator Pq,^^ has to be 
perturbed by an additional term of the form 

h^X^^ = h^e^^/^Xe-^^/^ = h^X - hf o ifXip. 

By (14. 3 p and the estimate ||/i^X^^t>|| < Cs/iHf ^, the inequahty (14.50 
may easily be perturbed into 

2(1 + C3)||Po,^.^ + h'X^Ml' > ^(1 - Cee)\\v\\' 
if h is small enough. Taking e small enough, we obtain 

(4.6) C,\\Po,^^v + h^X^M\'>^M' 

which implies with the choice of / = s + hs^ /2e that 

2 2 

Ce\\e'^eh\/\ + X)u\Y > h^\\e^ehu\f 
therefore we obtain the desired estimate since 

Remark 4.2. The use of Garding's inequality could give a stronger 
Carleman estimate, as in [26]. The present proof makes it possible to 
include boundary terms, which is useful in the study of inverse problems 
with partial data (see [H], [7] and [14]). 

In order to prove suitable solvability results, we need to shift the 
indices of the Sobolev spaces in the Carleman estimate by using pseu- 
dodifferential calculus. Recall that we assume that (M, g) is embedded 
in a compact manifold {N, g) without boundary, and that (fis a limiting 
Carleman weight near {U, g) where U is open in N and M U. The 
Laplace-Beltrami operator —A on N, with domain C°°(A^) C L'^{N), 
is essentially self-adjoint with spectrum in [0, oo). By the spectral the- 
orem we may define for s G R the semiclassical Bessel potentials 
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One has J^J* = J^^*, and commutes with any function of —A. 
Define for s G R the semiclassical Sobolev spaces via 

so H^^^{N) is the completion of C°°{N) in this norm. It is easy to see 
that the dual of H^^^{N) may be isometrically identified with H~J{N). 

It is a basic fact that J* is a semiclassical pseudodifferential operator 
of order s in iV (see [31] and [6]). This implies pseudolocal estimates: if 
ip,X ^ Co°{N) with X = 1 iisar supp ip, and if s,a, (3 E H and K eN, 
then 

(4.7) \\{1-x)J'Mh^jn) < CKh^'M^.^^^^y 
We will also use commutator estimates in the form 

(4.8) J>||i2(iV) < Ch\\u\\H^jN) 

whenever A is a first order semiclassical differential operator in A^. 

Lemma 4.3. Under the above assumptions on M, N, U and under the 
assumptions of Theorem 4jJ_, given s G R there are two constants Cs > 
and < kg < 1 such that for all functions u G C^{M°) and all 
< h < hs one has the inequality 

(4.9) l|e^M||//^=+i(Ar) < Csh\\e^{Ag + X + q)u\\Hs(N)- 
Proof. We consider the conjugated operator 

P^, = e^h\Ag + X + q)e-^ 

where (p^ is the weight defined in the proof of Theorem 14.11 Let x ^ 
C^{U) with X = 1 near M. Then the estimates (14. 6p . (14. 7p imply 

h\\u\\jjs+i < h\\xJ'u\\Hi + h\\{l - x)J'u\\h\ 

scl scl scl 

< CiVe\\P^,{xJ'u)\\L2 + Cih^\\u\\jjs+i. 

scl 

By the estimate || [P<o^, x]-^'*^||l2 < C2h'^\\u\\rTS+i, and by absorbing the 

scl 

error terms in the left hand-side if h is small enough, we obtain 

scl 

(4.10) h\\u\\j,s+. < CiV^\\rP^M\L^ + Civ^||x[P^., J>\\l^- 

scl 

Note that in estimating the last term one may extend if smoothly 
outside U if desired. Since 

P^^ = -h^A - iVv^eT + 2(V<^„ W) + HA'p, + h^X^^ + h^q 

and since [— /i^A, J^] = 0, the commutator estimates (14. 8 p imply that 
for ^ £ ^ 1 the last term in (14.101) may be absorbed in the left hand 

side. This finishes the proof since 
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Proposition 4.4. If A is a smooth l-form and q is a smooth function 
on M, there exists a constant < Hq < 1 such that for any func- 
tion f G L?'{M) there exists a solution u G H^{M) to the equation 
e'^/^Cg^A,q€-~'''^^u = f satisfying 

ll^llffi^iCAf) < Ch\\f\\L2(M)- 

Proof. We consider the conjugated operator 

Let / G L2(M), we consider the subspace E = P*{C^{N)) of H-l{N) 
and the hnear form defined on E by 

L{P;v) = {f\v)L^M), V G Co°°(M°). 

By Lemma 14.91 apphed to P* , this form is well defined and 

\LiP;v)\<\\f\\LHM)h-'\\P;v\\j,-i^j,y 

By the Hahn-Banach theorem, there is an extension L of L which is a 
bounded functional on H~l{N) with norm less than ||/||2,2. Since 
the dual of H~l{N) is H^^^{N), there exists a function u G H^cii^) 
such that L{v) = {u\v) and /iHmUji^i ^(^v) ^ ll/iU2(M)- Then u = u\m is 
the desired solution, since for all v G C^{M°) 

{p^u\v) = iu\p;v) = L{p;v) = L{p;v) = if\v). 

This completes the proof. □ 

5. Complex geometrical optics 

Let if he a limiting Carleman weight in an admissible manifold 
{M,g). We will construct solutions to Cg^gU = in M of the form 

(5.1) « = e-^(^+*'^)(a + ro). 

Here the real valued phase ip and complex amplitude a are obtained 
from a WKB construction, and the function tq will be a correction term 
which is small when h is small. 

We write p = if + iip ioi the complex phase. It will be convenient to 
extend the notations ( ■ , ■ ) and | ■ P to complex tangent vectors by 

(C, V) = (Re C, Re ri) - (Im C, Im rf) + i((Re C, Im ?]) + (Im C, Re 77)), 

icr = (c,c). 
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We make similar extensions for the inner product of cotangent vectors. 
With this notation, the conjugated operator Pp = e^^^h? Cg^qC'^/'^ has 
the expression 

Pp = -iVpP + h{2{yp) + Ap) + h'^Cg^q. 

Then (15. ip will be a solution of Cg^qU = provided that Pp{a + rQ) = 0. 
Following the WKB method, this results in the three equations 

(5.2) |VpP = 0, 

(5.3) (2(Vp) + Ap)a = 0, 

(5.4) Ppro = -h^Cg^qa. 

These equations will be solved in special coordinates in the admissible 
manifold {M,g). We know that {M,g) is conformally embedded in 
R X (Mo,5'o) ioT some compact simple (n — l)-dimensional (Mo,(?o)- 
Assume, after replacing Mq with a slightly larger simple manifold if 
necessary, that for some simple {D^qq) d (intMo,5'o) one has 

(5.5) (M, ^) (E (R X int D, ^) d (R x int Mq, g). 

Here R x Mq is covered by a global coordinate chart in which g has 
the form 



(5.6) g{x) = c{x) 



1 

go{x') 



where c > and go is simple. The limiting Carleman weight will be 

Lp{x) = Xi. 

5.1. The eikonal equation. Since ip was given, the eikonal equation 
(15.21) for the complex phase becomes a pair of equations for ip, 

|VV'P = |V<^P, (V<^,VV>) = 0. 

One has V9(x) = xi and the metric is of the form (15. 6p . so Vy^ = 
and iVv^l = ^- The eikonal equation now reads 

|V^| = -, d,,iP = 0. 
c 

Under the given assumptions on (M, g), there is an explicit construction 
for ijj. Let u; G -D be a point such that {xi,uj) ^ M for all xi. Denote 
points of M by X = {xi, r, 9) where (r, 9) are polar normal coordinates 
in [D^go) with center uj. That is, x' = expl^{r9) where r > and 
9 G S*"""^. In these coordinates (which depend on the choice of u) the 
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metric has the form 

/ 1 

g{x^,r,e) = c{x^,r,e) 1 

\ m{r,e) 

where m is a smooth positive definite matrix. 

To solve the eikonal equation it is enough to take ip{x) = ipuj{x) = r. 
With this choice one has p = Xi + ir and Vp = where 

d=-(— + i— 

2\dxi dr 

5.2. The transport equation. We now consider (15. 3p . In the coor- 
dinates {xi,r,9) one has 

and the transport equation reads 

49a + (^aiog^)a = 0. 

We choose a as the function 

a = |(^|-VV/2ao(xi,r)6(0) 
where da^ = and b(9) is smooth. 

5.3. Complex geometrical optics solutions. Finally, the equation 
( 15.41) may be written as 

This can be solved by using Lemma 14.41 We obtain tq satisfying 

\\r4HUM)<h. 
We record the properties of the solution just obtained. 

Proposition 5.1. Assume that {M,g) satisfies fl5.5p . fl5.6p and let 

q G C°°(M). Let u e D be such that {xi,uj) ^ M for all Xi. If {r,9) 
are polar normal coordinates in [D^qq) with center uj, then, whenever 
(9ao(xi,r) = and b{6) is smooth, the equation 

Cg^qU = in M 

has a solution of the form 

u = e-^("^+'^^) (|(7|-i/V/2ao(xi, r)6(^) + Tq) 
where ||ro||Hi (m) ^ h. 
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Next we consider the case where a magnetic field is present. The 
construction of complex geometric solutions is similar to the case with- 
out magnetic field, except that we have an additional factor in the 
amplitude. 

Proposition 5.2. Assume that {M,g) satisfies flS.Sp . fl5.6p . and let A 

be a smooth 1-form in M and q G C°°{M). Let u & D be such that 
{xi,uj) ^ M for all xi . If (r, 9) are polar normal coordinates in {D, go) 
with center u, then, whenever dao = and b is smooth, the equation 

Cg^A,qU = in M 

has a solution of the form 

^ = e-^("^+*")(|(7|"^/V/'e'*ao(xi,r)6(e) +ro) 
where $ satisfies 

(5.7) d^ + l{Ai + iAr) = 0, 

and ||ro||j^i^^(Af) < h. 

Proof. If Pp = eP/'^h'^ Cg^A,qe~^^^ I one computes 

Pp = -|VpP + h{2\/p + Ap + 2i{dp, A)) + h'^Cg^A,q- 

li p = Xi + ir then |Vpp = 0, and the transport equation for the 
amplitude a will be 

Ada + d(^\og ^-^^a + 2i{Ai+iAr)a = 0. 

A solution is given by a = \g\~^^'^c^^'^e'^'^ao{xi,r)b{6), where $ is a 
solution of fl5.7p and dao = 0. The equation for rg becomes 

and Lemma [4.41 finishes the proof. □ 

6. Uniqueness results 

Let (M, g) be an admissible manifold. In this section we will prove 
the global uniqueness results, Theorems [3] to [H As often in inverse 
boundary problems, the starting point is an integral identity which 
relates the boundary measurements to solutions inside the manifold. 

Lemma 6.1. Let Ai, A2 be smooth 1- forms in M and qi,q2 G C°°{M). 

U -^g,Ai,qi = -^g,A2,q2' thcn 

(6.1) / [i{Ai-A2,udv-vdu) + {\Ai\'^-\A2\^ + qi-q2)uv]dV = 0, 
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for any m, t> G H^{M) satisfying Cg^Ai,qiU = and Cg-A2,q2'^ = in M. 

Proof. One has A*^^ ^ = A^^^ ^ and 

{J^g,A,qf\h)dM = {dAu\dAv) + {qu\v), 

whenever Cg^A,qU = in M and u\sm = f, v\sm = h. These facts imply 
the identity 

((Ag,Ai,gi - A<,,A2,g2)/l^)aM 

= / [i{Ai - A2,udv - vdu) + {\Ai\^ -\A2\'^ + qi - q2)uv]dV 

J M 

where m, f are H^{M) solutions of Cg^Ai,qiU = and '^^,^2,92'^ = in 
M, which satisfy u\qm = f and v\qm = h. The result follows. □ 

6.1. Recovering a potential. Suppose that we have two potentials 
gi,g2 £ C°°(M) for which the corresponding DN maps are equal. Ac- 
cording to Lemma 16.11 one has 

(6.2) / {qi-q2)uiU2dV = 

J M 

for any Uj G H^{M) which satisfy Cg^q.Uj = in M. We use Proposition 
15.11 and choose solutions of the form 

u, = e-^(^^+*'-)(|(7|-i/V/2e^"(^^+*'^)6(^) +ri), 

where A is a real number and ~ ^- Inserting these solutions 

in (16.21) and letting h ^ shows that 

^iXix^+ir)^^^ _ q^)c{xi,r,e)b{e) drdOdxi = 0, 

with = {(r, 6*) ; {xi,r,0) G M}. We can extend qi — q2 smoothly 
by zero since qi = q2 up to infinite order on dM by Theorem 18. 4[ 
and we may then assume that the integral is over H x D. Taking the 
Xi-integral inside and varying b gives 

J e"^" J e^^^n^i - q2)cixi,r, 6) dxi^ dr = 0, for all 9. 

We denote the expression in parentheses by /(r, 9), and obtain that 




/(7(r))exp 



— / Xds 



dr = 
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for all D-geodesics 7 issued from the point uj. Varying uj and using the 
injectivity of the attenuated geodesic ray transform given in Theorem 
17.11 with constant attenuation —A, we obtain 



which holds for small enough A. Since [(l\—(l2)c is compactly supported 
in xi, its Fourier transform is analytic, and we obtain that q\ = q^- This 
proves Theorem [31 

6.2. Recovering a magnetic field. Next we show that also a mag- 
netic field can be recovered from the DN map. We will need the follow- 
ing standard reduction of the problem to a simply connected domain. 

Proposition 6.2. Let Mq, M he compact manifolds with boundary with 
Mq (e M, and let {g,Aj,qj) be smooth coefficients in M such that is 
not a Dirichlet eigenvalue of Cg^Aj,qj in Mq (j = 1,2). Suppose that 



Proof. Since the DN maps coincide in Mq, one has by Lemma [6.11 



/ [i{Ai- A2,UQdvQ-VQduQ) + {\Ai\'^ -\A2\'^ + qi- q2)uQVQdV] = 0, 

J Mo 

for any H^{Mq) solutions of Cg^Ai,qiUo = and £3,-^2,92 "^o = in Mq. 
If u and V are as above, then the restrictions to Mq solve the corre- 
sponding equations in Mq, and we obtain (16. ip since the coefficients 
coincide outside Mq. □ 

Proof of Theorem By boundary determination. Theorem 18.41 and 
after a gauge transformation, we may extend Aj and qj smoothly so 
that Ai = A2 and qi = q2 outside M. We are now in the setting of 
Proposition 16. 2[ Therefore, replacing M by a larger manifold inside 
R X Mq if necessary, we may assume that M is convex and Aj and qj 
are compactly supported in M, and also that the identity in Lemma 
16.11 holds whenever u,v E H^{M) are solutions of Cg^Ai,qiU = and 
Cg-A2,q2'^ = in M. 

Use Proposition 15.21 to choose solutions of the form 




^^{qi - q2)c{xi,r, 6) dxi = 0, for all r, 9, 



Ai = A2, qi = q2 in M\ Mq. 

If the DN maps Ag^Aj,qj in Mq coincide, then the integral identity (16.11) 
is valid for any H^{M) solutions of Cg^A^^q^u = and Cg^_A2,q2'^ = in 



M. 



u = 
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where p = xi + ir, dao = 0, and $ = $i + $2 is a solution of the 
equation 

(6.3) 9$ + = in M. 

Here A = Ai — A2, and Ai and Ar are the components in the xi and r 
coordinates. Inserting these solutions in (16.11) . multiplying both sides 
by h, and letting h ^ implies that 

lim / {A,dp)uvdV = 0. 



M 



Writing the integral in local coordinates gives 

(Ai +iAr)e''^ao{xi,r)b{9) dxi dr d9 = 0. 
Varying b leads to 

{Ai + iAr)e''^ao{xi,r) dp A dp = 0, for all 9, 



where fie = {(xi,r) G ; {xi,r,9) G M} is identified with a domain 
in the complex plane, with complex variable p. Integrating by parts 
and using (16.31) gives 

(6.4) [ e^*ao dp = 0. 

Jang 

The arguments in ^ Section 5], see also Section 7], then imply 
that e**|gng = F\g^g where F G C{Q0) is a nonvanishing holomorphic 
function, and F = where G G C{VL0) is holomorphic and G — 2$ is 
constant on dVLa. We choose 



ao 



where A is a real number, and then (16. 3p . (16. 4p . and integration by 
parts imply 

! (ii + ii,)e^^(^i+^") dp A dp = 0, for all ^. 

We define 

/(x')= / e'^^Mi(xi,x')dxi, 



a{x') = J2{ e'^'''^j{xux')dxiyx^. 

j=2 
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The identity above may be written as 

j e"^" [f{l{r)) + ia(7(r))] dr = 0, for all 9. 

The r- integrals are integrals over geodesies j in D. By varying the 
point uj in Proposition 15.21 on dD and using Theorem 17.11 we see for A 
small that / = —Xp and a = —idp where p G C°°{D) and p\dD = 0. 
The definition of a and analyticity of the Fourier transform imply that 

dkAj - djAk = 0, j, A; = 2, . . . , n. 

Also 

j e'^'^'idjAi - diAj){xi, x') dxi = djf + iXa^ = 0, 

showing that dA = in M and the magnetic fields coincide. 

Since M is convex and A\gM = 0, one has A = dp where p G C°°(M) 
can be chosen so that p\dM = 0. By a gauge transformation, we may 
assume that Ai = A2, and this 1-form will be denoted by A. The 
integral identity (16. ip becomes 

(6.5) / iqi-q2)uvdV = 0, 

for any solutions of Cg^A,qiU = and Cg^^A,q2'v = 0. Use Proposition 
15.21 and choose solutions 

where p = Xi + ir, and $ satisfies 

d<^ + ^{Ai + iAr) = inM. 
Letting /i — > in (16. 5p gives 

Proceeding as in the proof of Theorem [3] shows that qi = q2- D 

6.3. Recovering a conformal factor. The results on the Schrodinger 
inverse problem can be used to recover a conformal factor from the DN 
map. Recall that we use the notation = A^ when the potential q 
is zero. 

Proof of Theorem\^ It is enough to show that if {M,g) is admissible 
and c is smooth and positive, and if A^g = Ag, then c = 1. We have 
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c\dM = 1 and duc\sM = by Proposition 18. 11 and then the assumption 
and Proposition 18.21 imply 

where q = — Ag(c^)/c'^. We conclude from Theorem [3] that g = 0, 
so Ag(c^) = in M. Since = 1 on 5M, uniqueness of solutions 
for the Dirichlet problem shows that c = 1. □ 



7. Attenuated ray transform 



For the uniqueness results in inverse problems, we have used that 
certain geodesic ray transforms are injective. If (M, g) is a compact 
manifold with smooth boundary, geodesies can be parametrized by 
points on the unit sphere bundle SM = {{x,C,) G TM; |^| = 1}. For 
(x, ^) e SM let 7x,^(^) be the geodesic with 7(0) = x and 7(0) = ^. We 
assume that {M,g) is nontrapping, which means that the time t{x,^) 
when 7^.^^ exits M is always finite. 

Given a smooth real function a on M, the attenuated geodesic ray 
transform of a function / is given by 



/(7x,^(t))exp / a{-f^^^{s)) ds 



dt 



for (x,^) e d+SM. Here we use the sets of inward and outward pointing 
unit vectors 

d±S{M) = {{x,0 eSM;xedM, ±{^,iy{x)) < 0}, 

and u is the outer unit normal vector to dM. 

We will also need to integrate 1-forms over geodesies. Let / be a 
smooth function and a = ttj dx^ a smooth 1-form on M, and consider 

(7.1) F{x,0 = f{x)+a,{x)C 

for (x, ^) G SM. The attenuated geodesic ray transform of F is 



PF{x,0 



F{lx4{t),'yx,i{t)) exp / 0(7^,5(5)) ds 



dt. 



This transform always has a kernel: if p is a smooth function on M 
with p\dM = 0, then a direct computation shows that 

r{ap + dp{O)=0. 



The main result in this section states that for simple manifolds and 
small attenuation, this is the only obstruction to injectivity. 
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Theorem 7.1. Let {M,g) be a compact simple manifold with smooth 
boundary. There exists e > such that the following assertion holds 
for any smooth real function a on M satisfying \a\ < e: If F is given 
by dZH) and if 

for all {x,^) G d+SM, then F = ap + dp{^) for some smooth function 
p on M with p\dM = 0. 

Note that if a = 0, this shows that any function / whose attenuated 
ray transform vanishes must be identically zero. If / = and a ^ 
everywhere, then any 1-form whose attenuated ray transform vanishes 
must be identically zero. Injectivity of the geodesic ray transform for 
functions and 1-forms on simple manifolds in the case a = is well 
known [1], [20], [27]. The injectivity for functions and small a is proved 
in [27], [2H] under conditions which involve a modified Jacobi equation 
or the size and curvature of M. We give a proof which works on simple 
manifolds. 

We remark that the notation in this section is somewhat different 
from the other sections. For instance, we will denote by V the covariant 
derivative and more generally the horizontal derivative. The notation 
will be explained in more detail below. 

7.1. Preliminaries. The proof will be based on energy estimates and 
a Pestov identity, which is the standard approach to such problems. 
First we need to recall the definition of horizontal (or semibasic) ten- 
sor fields on TM. These are {p, q) tensor fields on TM which have 
coordinate representations 

• • • • d d 

with respect to coordinates (x, ^) on TM associated to charts x of M. 
The components transform in the same way as tensors on M under 
changes of charts. Tensor fields on M may be considered as ^-constant 
horizontal tensor fields on TM. See [27] for an invariant definition and 
other details on horizontal tensor fields. 

For our purposes, it is sufficient to know that smooth functions on 
TM are horizontal tensors of degree (0, 0), and that the horizontal and 
vertical derivatives defined by 
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are invariantly defined operators which map horizontal {p, q) tensors to 
horizontal {p + l,q) tensors. Here 

f) " " 

Thus V acts in the same way as the usual covariant derivative. 

Below, we will work with smooth functions and tensors on SM. Let 
p : TM \ {0} — >• SM be the map (x,^) h- >• (a;, ^/|^|). The horizontal 
and vertical derivatives on SM are defined by 

h V 

Vu^V{uop)\sM, du^V{uop)\sM- 

h V 

Also ViU = Vj(ti o p)\sM and diU = Vi{u o p)\sM- H u is horizontal 
tensor field on SM, then and 9ii are also horizontal tensor fields. 
For a smooth function on SM, these derivatives are given by 

We mention the following formulas 

[Vj, dj] = 0, [di, dj] = .^jaj- - ^jdi, 

where R is the curvature tensor and m is a scalar function. We write 

d^u — g^^djU, {du, dv) — d^udiV, \du\^ — d^udiU. 

Let H be the geodesic vector field on SM that generates geodesic 
flow. In local coordinates 

Hu{x,^) = c^iu{x,^/m - r\,{x)eediu, e e 

where -u is a smooth function on SM. We may apply the operator H 
to horizontal tensor fields on SM by defining Hu = ^^ViU. 

7.2. Ray transform of functions. We now consider the boundary 
value problem for transport equation: 

Hu + au^ - f, u\d_s{M)^0, 
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where the absorption a and / are smooth functions on {M,g). The 
solution is given by 

dt. 



ui{.x,i)= I /(7^,5(t))exp / a(7^,5(s)) c/s 
Jo L Jo 

The trace 

^"f = ui\d+s{M) 

is the attenuated geodesic ray transform of the function /. It is natural 
to define T\g_sM = 0, and then ul indeed vanishes on d-SM. 

We will prove that / is uniquely determined by 1°"/ under the follow- 
ing assumption. Let 7 = 7^,^, (x, ^) G d+SM, be an arbitrary geodesic, 
and consider the quadratic form 



(7.2) 



E^{X)= [ {\DX\'^ - {R^X,X) - a^\X\'^){t)dt 
Jo 



where r = r(x,^), X{t) is a vector field on 7 belonging to the space 

i/o'(7) = {X e H\[0, r]; r(7)) ; X(0) = X(r) = 0}, 

D is the covariant derivative along 7, and R^X = R{X, 7)7. We assume 
for all geodesies 7 the positive definiteness of this quadratic form, 



(7.3) 



E^iX) > whenever X e H^{-f), 
E^{X) =0 iff X = 0. 



If (17. 3p holds we say that any geodesic has no conjugate points with 
respect to (17.21) . If a = we obtain the usual index form E^ = E^. 
Then clearly there are no conjugate points in the usual sense if there 
are none with respect to (17. 2p . 

Proposition 7.2. Let {M,g) be compact and dM strictly convex. As- 
sume that any geodesic has no conjugate points with respect to i \7.2L 
Then any smooth function f on the manifold (M, g) is uniquely deter- 
mined by its attenuated geodesic X-ray transform. 

Proof. Let /"/ = 0. We will assume that m = m{ is smooth on SM 
(otherwise one can work in a slightly smaller manifold than M, and 
pass to the limit using the smoothness properties of r as in [27]). The 
function u satisfies 

dHu + adu = 

and therefore 



(7.4) 



42 LIMITING CARLEMAN WEIGHTS 

Using the formulas for V and 9, it is not difficult to check the following 
identity (valid for any u G C°°(5'M)): 

(7.5) \dHu\^ = \Hdu\^ + 6V + OW - R{du, ^, ^, du), 

where S and 6 are the vertical and horizontal divergences, 

6X = ViX\ ex = diX\ 

and V and W are defined by 

V' = {du, Vu)C - {Hu)d'u, W' = {Hu)V'u. 

From (17. 4p and ( 17.5^ we obtain 

\Hduf - R{du, ^, ^, du) - a'^\du\'^ + SV + eW = 0. 

Now integrate this equality over the manifold SM. Before this we recall 
the integration formulas (see |27j): 



V d{SM) = dM I vdS, 

SM J M 



XI 



eXdS, = {n-l) / {X,i)dS^ 

Sx ^ Sx 

V / udSr = / VudSr, 



X 

Sx " Sx 



where w is a scalar, X is a horizontal vector field and u is horizontal 
tensor field, 

/ 6XdM= I {X,u)d{dM), 

Jm JdM 

where X is a vector field on (M, g). In these formulas the volume forms 
of corresponding manifolds are naturally defined using the metric g. 
After integration we have 

(7.6) f {\Hdu\'^ - R{du,^,^,du)-a^\du\'^)d{SM) 

JsM 

+ {n-l) I \Hu\'^d{SM) = 0. 

JsM 

We used the fact that (V, z/) vanishes on d{SM) since u\q(^sm) = 0. 
We next show that (17.31) implies 

E%Y)= [ {\HY\'~R{Y,^,^,Y)-a'\Y\')d{SM)>0, 

J SM 

for any horizontal vector field Y G C^{SM]TM) with Y\q(^sm) = 0, 
and that the equality holds iff F = 0. Santalo's formula (see [27j ) 
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states 

/ vd{SM) = - [ '^\{^,,^{t),^^,i:{t)){^,iy) dtd{d{SM)), 

J SM Jd+S{M) Jo 

for V e C°°(^M). Let 

X(x,e,t) =r(7,,^(t),7,,^(t)), 
which imphes HY{'yx,iit),'jx,^(t)) = DX{x,C,,t). Then we have 

E^iY) = - [ El {X){x,0{^,'y{x))d{d{SM)) > 0. 

Jd+S{M) 

Equahty holds iff F = 0. We have from (17.61) that Hu = 0, which 
imphes u = and / = 0. □ 

7.3. Ray transform of 1-forms. Let / be a smooth function and 
a = ai{x)dx'^ a smooth 1-form in M, and consider the boundary value 
problem 

Hu + au = -F, u\o_sM = 0, 
where a is a smooth function on M and F is as in (17.11) . The solution 
u = is given by 

rT{x,i) r-t 



dt. 



and the trace 

PF = U^ld^SM 

is the attenuated geodesic X-ray transform of F. 



Proposition 7.3. Let {M,g) be compact and dM strictly convex, and 
suppose that any geodesic has no conjugate points with respect to i\7.2\\ . 
If I'^F = 0, then F = ap + dp{^) for some smooth function p on M 
which vanishes on dM. 

Proof. We follow the proof of Proposition 17.21 Again assume that u is 
smooth in SM. Then u satisfies 

dHu + adu = —dF, 

and 

\dHu\^ = a^\du\^ + 2a{du, OF) + \dF\\ 
The identity (17.51) then implies 

\Hdu\'^-R{du,^,^,du)-a'^\du\'^ + 6V + eW-2a{du,dF)-\dF\'^ = 0. 
Integrating this over SM, it follows that 

(7.7) E''{du)+ [ [{n-l)\Hu\^ -2a{du,dF) -\dF\^]d{SM) = 0. 

J SM 
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Using the integration formula 

/ divdS^ = (n-l) / v^idS^ 

and the identities 

\Huf = {auf + 2auF + 

we easily obtain from (17.71) that 

E^idu) + (n - 1) / {au + ff d{SM) = 0. 

JsM 

The assumption (17. 3p and Santalo's formula imply E"-{du) > 0, so 
from the last equality we obtain E°'{du) = and au + f = 0. Now 
E°'{du) = implies du = 0, so u = u{x) and 

{du,0 + au = -f - {a,0, xeM,^eS^. 
The claim follows. □ 

It remains to show that simple manifolds satisfy the condition of 
Proposition 17.31 



Proof of Theorem 7.1. Let {M,g) be a compact simple manifold with 
smooth boundary. We need to show that there is e > such that 

(7.8) E,^^^{X)>e r'^^'^\x\'dt, 



for all (x,0 e d+SM and for all X e H^ilx,^)- If UM holds and 
\a\ < ^/e then any geodesic on {M,g) has no conjugate points with 
respect to (17. 2p . and Proposition 17.31 implies the desired result. 

Let first (x,^) e d^SM, and consider the unbounded operator on 
L^{'lx,(,) with domain fl HQi^j^^^), given by 

L^^ ^ : X ^ -D'X - i?(X,7,,^)7,,5. 

This operator is self-adjoint and has discrete spectrum, which lies in 
(0, oo) since the corresponding quadratic form E^^ ^ is positive definite. 
Therefore 

E^^^iX) > Ai(a;,0 [ ^ '^\x\'dt, X G H^oilx^^ 
Jo 

where Ai(a:, ^) > is the smallest eigenvalue. The coefficients of L 
depend smoothly on (x, ^) and r is smooth and positive in d^SM, so 
it is not hard to see that (17.81) holds in a neighborhood of any fixed 
point in d+SM. 
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For tangential directions we use the Poincare inequality 

f \uit)\^dt<^ [ \u{t)\^dt, ueH^{[0,L]), 
Jo ^ Jo 

where the constant L^/tt^ is optimal [23j. If (x, ^) G S{dM) and 6 > 0, 
then by continuity of r there is a neighborhood U of (x, ^) in d+SM 
such that r(?/, r/) < 5 in that neighborhood. Choosing 5 small enough, 
the Poincare inequality implies 



Jo 



whenever (y,?7) G t/ and X G ifo(7y^^). This shows (17.81) near any 
point of S{dM). It follows that for some e > 0, (17.81) holds on the 
compact set d+SM. □ 

8. Boundary determination 

To deal with the inverse problems we are interested in, we need to 
use the fact that the DN map determines the Taylor expansions at 
the boundary of the different quantities involved. In the case of the 
Laplace-Beltrami operator, the relevant result is proved in [18] and is 
as follows. 

Proposition 8.1. Let {M,gi) and {M,g2) be compact manifolds with 
smooth boundary, with dimension n > 3. If Ag-^ = Ag^, then the Taylor 
series of gi and g2 in boundary normal coordinates are equal at each 
point on the boundary. 

In this section we extend the previous result to the case where electric 
and magnetic potentials are present. First we need to consider the 
gauge invariance of the DN map. 

Proposition 8.2. Let {M,g) be a compact manifold with boundary, 
and let A be a smooth 1-form and q a smooth function on M . If c and 
i/j are smooth functions such that 

c> 0, c\aM = 1, duc\aM = 0, ipldM = 0, 



then we have 

7 n — z ^ / n — 2 \ 

where gc = c * (c 4 j . 



g,A,q -^^c~^g,A+dtl>,c{q—gc) 

n—2 ^ , n—2 s 
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Proof. It follows from a direct computation that 

C 4 Lg^A,q{C 4 U) = Lc-lg,A,c{q-qc)'^, 
e ^^Cg^A,q{e^^u) — Cg^A+dil>,qU- 

Let / e C°°{dM), and let u be the solution of Cg^A,qU = in M which 
satisfies u\qm — f ■ If 



and if M = c * e ''^m, we have -m = in M and mI^m = /• Then 



by using the boundary values of c and and the fact that i/g — Ug. □ 

Remcirk 8.3. The conformal gauge invariance is related to the fact 
that the conformal Laplace-Beltrami operator Ag on {M,g), defined 

by 

4(n — 1) 

transforms under a conformal change of metrics by 

We use the notation /i ~ /2 to denote that /i and /2 have the same 
Taylor series. Our main boundary determination result is as follows. 

Theorem 8.4. Let {M,gi) and {M,g2) be compact manifolds with 
boundary, of dimension n > 3, and let Ai,A2 be two smooth 1-forms 

and qi,q2 two smooth functions in M. If Ag-^^Ai,qi = ^92,^2,92 ^'^^ ^/ 
p e dM, then there exist smooth positive functions Cj with 

CjldM^i and duCjloM^O, 

and smooth functions ipj with ipjlaM — 0, such that the gauge trans- 
formed coefficients 

{9j, Aj, qj) = {cfgj, Aj + dipj, Cj{qj - dfAg.dj)) 

71-2 

with dj — Cj * , satisfy in boundary normal coordinates at p 

gi ~ g2, Ai ~ A2, qi ~ q2- 

Furthermore, if gi ~ g2 in boundary normal coordinates on all of dM, 
then Ai ~ A2 and qi ~ q2 on dM. 
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Boundary normal coordinates (x', x„) at a boundary point p are such 
that x' parametrizes dM near p and x„ is the distance to the boundary 
along unit speed geodesies normal to dM. See [18] for more details. 

We prove Theorem 18.41 by showing that Ag^A,q is a pseudodifferential 
operator whose symbol determines the boundary values of the coeffi- 
cients. This method was used in [33] for the conductivity equation, 
in [18] for the Laplace-Beltrami operator, and in [22] for the magnetic 
Schrodinger operator with Euclidean metric. We follow the proof in 

m- 

The first issue to consider is the gauge invariance in the operator 
^g,A,q- This will be dealt with by normalizing the coefficients {g, A, q) 
in a way which is suitable for boundary determination. Fix a point 
p G dM and boundary normal coordinates [x', Xn) near p. In these 
coordinates dM corresponds to {x„ = 0} and 

9 = 9ai3 dx"' (g) dx'^ + dx^ ® dx"" . 

We use the convention that Greek indices run from 1 to n — 1 and 
Roman indices from 1 to n. 

Lemma 8.5. Let {g, A, q) be smooth coefficients in M and p G dM. 

There exist a positive smooth function c with c\qm = 1 and d^c\gM = 0, 
and a smooth function ip with ip\dM = 0, such that in boundary normal 
coordinates near p the quantities g = c~^g and A = A + dip satisfy 

(8.1) diAn{x',0)=0, j>0, 

(8.2) dii~go.pdnr''){x',0)=0, j>l. 

Proof. One can find a smooth function ip near p with tp{x', 0) = and 
dl'^^ip^x' ,0) = —dlAn{x',0) for j > 0. Extending this in a suitable 
way, one obtains ip G C°°(M) with iP\qm = such that A = A + dip 
will satisfy flO) . 

Further, we construct a smooth function c near p which satisfies 
c(x', 0) = 1, dnc{x', 0) = 0, and 

^(logdet(c^7"^))(x',0) = 0, j>2. 

In fact, one may take c = where /i(x', 0) = dnn{x', 0) = and 
(9^/i(x',0) = -^a^(log det(^"^))(x',0) for j > 2. There is an ex- 
tension of c to a positive function c G C°°(M) with c\qm = 1 and 
duc\dM = 0. Since 

a„(log detir^)) = -g^pdnr^ 



one also has the second condition (18.21) . 



□ 
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Replacing coefficients {g, A, q) by the gauge equivalent coefficients 
{c~^g, A + dip, c{q — d~^Agd)) does not affect the DN map. Therefore, 
below we will assume that (18. ip and (18.21) are valid. Note that boundary 
normal coordinates for g are also boundary normal coordinates for any 
conformal multiple of g, if the conformal factor is 1 on dM. 

The next aim is to write Ag^^^g as a pseudodifferential operator and 
to compute the symbol in a small neighborhood of p. Here we use the 
usual (not semiclassical) pseudodifferential calculus and left quantiza- 
tion, so that a symbol r(x,^) in T*R" corresponds to the operator 

Rf{x) = (27r)-" / / e'^'-y>^r{x,Of{y) dyd^. 

We denote by p ~ J2Pj asymptotic sum of symbols, see [10] for 
these basic facts. 

Lemma 8.6. Ag ^ g is a pseudodifferential operator of order 1 on dM. 
Its full symbol (in left quantization) in boundary normal coordinates 
near p is —h ^ — X]j<i where bj are given in (18. 5p - (18. 8p . 

Proof. In the x coordinates, one has 

Cg^A,, = -Ag + 2g^''A,Dk + G, 

where G = \g\-^/^Dj{\g\^/'^g^''Ak) + g^^AjAk + q and 

\g\ = det{gjk) = detig^fj). 

From (18. ip we know that d^An{x', 0) = for all multi- indices K G N". 

One would like to have An = also inside M. To achieve this, we 
introduce as in [22] the conjugated operator 

M = e-'^Cg,A,,e'\ 

where h{x) = - J^" A^ix' , s) ds. Note that d^h{x\0) = 0. Writing 
Aj = Aj + djh, we obtain An = and 

M = -Ag + 2g'"'A^Dp + G, 

where G = \g\-^/^Da{\g\^^^g"^Ap) + g°''^AaAp + q. We then have 

M = Dl + iE{x)Dn + Q2{x, D,,) + Qi{x, D,,) + 2g''^A^Dp + G, 

with E, Qi, and Q2 given by 

E{x) = ^g^pdng''^, 
Q2{x,D,.)=g''^D^Dp, 
Qi(x, D,,) = -r{lg^Pd^{\og \g\) + do.g''^)Dp. 
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As in [18], one would like to have a factorization 



(8.3) M = {Dn + iE{x) - iB{x, D,,)){Dn + iB{x, D^,)) 

modulo a smoothing operator, 

where S is a pseudodifferential operator of order 1 with symbol ^'). 
Using left quantization for symbols, (18. 3p implies on the level of symbols 
that 

(8.4) dnh-E{x)h+Y^ = g2 + gi + 2(?"^iae/3+G mod 5-°°. 

K 

Let h ~ X]j<i "where hj{x, is homogeneous of order j in Insert- 
ing this in (18. 4p and collecting terms with the same order of homogene- 
ity, one obtains hj as follows: 

(8.5) h = -v^, 

(8.6) 6o = ^{-dnhi + Eh - V^^i ■ D^,hi + q^ + 2g^''A^ip), 

1 ^ d^biD5bk ~ 

(8.7) h^, = —{-dr.h, + Eh,- J2 ' ^, 

^ 0<j,fc<l 
j+k=\K\ 

(8.8) b^^, = —{-dnbm + Ebm- J2 K\ ^ i^<-^)- 

^ m<j,A;<l 
j+k—\K\=m 

With these choices, b E and (18. 3p is valid. By the argument in fl8\ 
Proposition 1.2], one has 

Ag,Aji^') = -B{x',0,D,,)f{x') + Rf{x') 
where i? is a smoothing operator. □ 

From the symbol of the DN map, one can recover the following in- 
formation on the coefficients. 

Lemma 8.7. The knowledge of ^g,A,q determines on {xn = 0} the 

quantities 

(8.9) dng"'^, d^Ao,, d^T^. 
Here K G N" is any multi-index and 

with = 9„(7"^ - {g^5dng^^)g''^ ■ 
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Proof. By Lemma ESI ^9,A,q determines &|x„=o and each bj\x„=o- From 
bi\x„=o one recovers g"'f^\x^=Q. If T{g°'^) denotes any linear combination 
of tangential derivatives of g""^, one has for x„ = 

bo = ^(-^n^i + Eh + 2g^PA^ip) + T{g^^) 
= -^-f^ + \E-^g-n^i, + T{g'^P) 

4g2 2 ^ 

= —k'^^u^up - g'-n^up + T((7°^), 

where u = C/\C\g- Evaluating at ±lj and varying u one recovers Aa 
and and consequently also dng"'^, for x„ = 0. 
Moving to 6-1, one has for x„ = 

6-1 = 7^{-dnbo + Ebo + G)+ Tig'^", dng''^ i«) 

where T denotes tangential derivatives of the given quantities. Thus, 
one recovers /"^ and dnA^ on {xn = 0}. By induction, we prove that 
for J > 1 

h-, = - ( - ^)'((5r^r^)c^au;/3 + 9'''{diA^)u:p) 

Indeed, this is true for j = 1, and assuming this for j one gets 

= -^dnb-3 + Tig'^^ dng'^^ r\ d^H^^ i„, . . . , diA^) 
= - ( - ^y\{dT')^.uj, + g-^{di^^A^)u^,) + r( ■ ). 

Thus one recovers 9^/°'^ and c?^Aq, on {x„ = 0} for all j > 0. The 
result follows since Aa = Aa when x„ = 0. □ 

We may now prove the main result on boundary determination. 

Proof of Theorem \8.4\ Let Ag A,q = Ag- Replacing both sets of co- 
efficients by gauge equivalent ones as discussed after Lemma 18.51 we 
may assume that (18.11) and (18. 2p are valid. Then Lemma 18.71 implies 
that the quantities (18. 9p with and without tildes coincide on {xn = 0}. 
We prove by induction that for j > 0, one has on a;„ = 

(8.10) dig = dig, 9^+V^ = a^+V^ gapdi^'g^^ = gapdt'r^. 
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We first note that 

(8.11) go.pdlg'^^ = -g^pdlr^ on x„ = 0. 

This follows from (18. 2p for g and g, since d^g"'^ = d^g"'^ on x„ = for 
j = 0, 1. Note also that (18.21) implies 

dik'^^ = di+'g"^ - {g,sdng'')V^ on x„ = 0, 
and therefore 

= lidi-^'g'^^ - {9,8dn9^')di^'9''^) + di{qg'^^) on x„ = 0. 

To prove (18.101) for j = 0, we use that 1°'^ = on x„ = 0. This 
implies, upon multiplying by gap and summing, that g = g on x„ = 0. 
Here we used (18.111) . Then = also implies d^g"^ = d'^g'^^ on 
= 0. The equality ga/sdnO"^ = Qapd^g""^ follows by using (18.21) . 

Assume now that ( ]8.10p holds for j < k. Moving to A; + 1, the 
equality 9^+^/°^ = d^'^ 1°^ on x„ = imphes upon multiplying by gap, 
summing, and using the induction hypothesis, that 

g.pdt^qg''^) = ~go.pdt\qr^) on x„, = 0. 

The induction hypothesis again gives d^^q = d^^q on x„ = 0, and 
gk+3ga/3 _ Qk+3ga/3 |;]^g^ follows by goiug back to the equality Q^"*"^/"^ = 
Qk+ij^ap _ g_ ^Y^Q lagt statement in (18.101) for j = A; + 1 is a 

consequence of (18. 2p . This ends the induction. 

The outcome of the above argument is that (7"^ ~ — A^, 

and g ~ g at p. This shows the first statement in Theorem 18. 4[ If ~ ^ 
at each p, then it is easy to obtain g ~ g at each p from 1°'^ ~ P^. 
Also, the function ip constructed locally in Lemma [831 can be obtained 
globally on dM by a suitable partition of unity. Therefore A ~ A on 
all of dM. □ 

Appendix A. Riemannian geometry 

In this appendix we include basic definitions and facts which are used 
throughout the text. For more details see [12j. We are using Einstein's 
summation convention: repeated indices in lower and upper position 
are summed. In the following (M, g) is a Riemannian manifold. When 
no confusion may occur, we will use the following standard notations 
for the inner product and the norm: 



{X,Y)=giX,Y), \X\ = y/g{X,X). 
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A.l. Connection and Hessian. The Riemannian metric g on M in- 
duces a natural isomorphism between the tangent and cotangent bun- 
dles given by 

T(M) ^ T*{M) 
(x,X) ^ (x,X^) 

where X^(Y) — {X, Y), and whose inverse is 

T*(M) ^ T(M) 

where is defined by ^{X) = {^'^,X). In local coordinates, if the 
metric is given by 

g — gjk dx-' (8) dx'', 

this reads 

X'^gjkX^dx'', C^^9'%d,,. 

In particular, the gradient field is defined by Vip = d(p^. The mu- 
sical isomorphisms allow to hft the metric to the cotangent bundle. 
The cotangent bundle is hence naturally endowed with the Riemann- 
ian metric g~^ given in local coordinates by 

g-^^g'^'dij^dik- 

It is natural to use (■, ■) and | ■ | to denote the inner product and the 
norm both on the tangent and cotangent bundles. 

We denote by D the Levi-Civita connection on (M, g). A connection 
is a bilinear map on the vector space of vector fields which satisfies the 
following conditions: 

(0 DfxY = fDxY, and Dx{fY) = {Lxf)Y + fDxY if / is a 

smooth function on M, 
(n) DxY-DyX=[X,Y]. 

Here Lx is the Lie derivative. On a Riemannian manifold, there is 
precisely one connection, called the Levi-Civita connection, which is 
consistent with the metric, i.e. which satisfies 

(m) Lx{Y, Z) = {DxY, Z) + {Y, DxZ). 

This connection is determined in local coordinates by 

where the Christoffel symbols F' j. are given by 
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Note that F^j^, = T\j. If X is a vector field on M, then the endomor- 
phism Dx has a unique extension as an endomorphism on the space of 
tensors satisfying the following requirements: 

(i) Dx is type preserving, 

(u) Dx{c{S)) = c{DxS) for any tensor 5* and any contraction c, 
(m) Dx{S ® T) = DxT ® 5 + T ® DxS for any tensors S, T. 

In particular, if / is a function we have Dx^p = dip{X) = Lx<p, and for 
1-forms the connection is given in local coordinates by 

The total derivative DS of S is the tensor DS{X, ■ ) = DxS{ ■ ). 

The Hessian of a smooth function is the symmetric (2, 0)-tensor 
D'^if = Ddip. The expression of the Hessian in local coordinates is 

dV = {dl^.,^ - r'^A^)dx^ ® dx'. 

The following identities will be useful: 

(A.l) DMX, Y) = h^^g{X, Y) = {DxV^, Y), 

(A.2) DMX,X) = ^^{^{t)) . 

Here 7 is the geodesic with 7(0) = X. 

A.2. Parallel and Killing fields. First we recall the following iden- 
tities for the Lie derivative: if / is a function and X a vector field 
then 

(A.3) Lfxg = fLxg + df®x' + x^® df, 

and if 5* is a (2, 0)-tensor then 

(A.4) iLxS)iY, Z) = Lx{S{Y, Z)) - S{[X, Y^Z) - S{Y, [X, Z\). 
Definition A.l. A vector field X in {M,g) is a Killing field if 

Lxg = 0. 

Note that flA.4D implies 

{Lxg){Y, Z) = {DyX, Z) + (y, DzX). 

Definition A.2. A vector field X is said to he parallel if its covariant 
derivative vanishes identically, that is DX = 0. 

The following characterization is used in the proof of Theorem [H 
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Lemma A. 3. Let {M,g) be a simply connected Riemannian manifold. 
A vector field on (M, g) is parallel if and only if it is both a gradient 
field and a Killing field. 

Proof. In a simply connected manifold, a vector field X is a gradient 
field if and only if the one form u — is closed. We have 

du;(Y,Z) = Ly{u;{Z)) - Lz{oj{Y)) -oj{[Y,Z]) 

= Ly{X, Z) - Lz{X, Y) - {X, [Y, Z]) 

= {DyX, Z) + {X, DyZ) - {DzX, Y) - {X, DzY) 

-{X,DyZ)^{X,DzY). 

Thus, X is a gradient field if and only if for all vector fields Y, Z 

{DyX,Z)-{DzX,Y)^0. 

On the other hand, X is a Killing field if and only if 

LxgiY, Z) = {DyX, Z) + {Y, DzX) = 

for all vector fields Y, Z. The result ensues. □ 

The next result states that the existence of a unit parallel vector 
field implies a local product structure on the manifold. 

Lemma A. 4. Let X be a unit parallel vector field in a manifold (M, g) . 
Near any point of M, there are local coordinates x such that X — d/dxi 
and the metric has the form 

Conversely, if such coordinates exist then X = d/dxi is unit parallel. 

Proof. Let X be unit parallel and let F be the distribution orthogonal 
to X. If V, are vector fields orthogonal to X then 

{[V,W],X) = {DvW - DwV,X) = V{W,X) -W{V,X) =0. 

Thus r is involutive, and by the Frobenius theorem there is a hypersur- 
face S (through any point of M) which is normal to X. Let x' i— >■ q{x') 
be local coordinates on S, and let {xi,x') i— > expg(^^,^{xiX{q{x')) be 
corresponding scmigcodcsic coordinates. In fact integral curves of X 
are geodesies (if ■j{t) = X{j{t)) then Da^j = 0), so X = c^i. Then 

H J) k >2 then digjk = {Dg^dj, d^) + {dj, Dg^dk) = since di is parallel. 
Therefore go = go{x'). 

The converse follows since Dg.di = by a direct computation. □ 
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Finally, wc need the definition of conformal Killing fields. 

Definition A. 5. A vector field X in (M, g) is called a conformal 
Killing field if 

Lxg = Agr. 

Note that if Lxg — \g.,hy taking traces one has A = | diw X. The 
notion of conformal Killing field is of course invariant under conformal 
change of the metric. 

A. 3. Submanifolds. Let (M, g) be a Riemannian manifold and let S 
be an embedded hypersurface in M . Fix a unit vector field which is 
normal to S. The second fundamental form of S is defined by 

i{XX) = {DxN,Y), 

where X and Y are vector fields tangent to S, and D is the Levi-Civita 
connection in (M, g) . The eigenvalues of the symmetric bilinear form 
I are the principal curvatures of S] their sign depends on the choice of 
normal. 

Definition A. 6. A point of a hypersurface is called umbilical if all the 
principal curvatures are equal at that point. A hypersurface is called 
totally umbilical if every point is umbilical. 

Definition A. 7. A hypersurface is called strictly convex if the second 

fundamental form is positive definite. 

A. 4. Curvature tensors. Next we consider curvature tensors on a 
Riemannian manifold {M,g). The Riemann curvature endomorphism 
is a (3, 1) -tensor on M, defined by 

R{X, Y)Z = DxDyZ - DyDxZ - Dyx,Y]Z, 

whenever A, F, Z are vector fields on M. By lowering indices, one 
obtains the Riemann curvature tensor which is the (4, 0)-tensor 

R{X, r, Z, W) = {R{X, Y)Z, W). 

In local coordinates, with coordinate vector fields da = d/dxa and with 
Da — Dq^ , the components of the curvature tensor are given by 

Rabcd = {{DaDh - DbDa)dc, da). 

By taking traces of the Riemann curvature tensor, we obtain the Ricci 
tensor which is a symmetric (2, 0)-tensor whose components are 

Rbc — g"''^ Rabcd- 

The scalar curvature is the function 

Seal = y^^i?6c. 
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A manifold (M, g) is said to be flat if the Riemann curvature tensor 
vanishes identically. Euclidean space is flat, and any flat manifold is 
locally isometric to a subset of Euclidean space. 

Finally, we wish to define the conformal curvature tensors. First con- 
sider the rho-tensor, which is a symmetric 2-tensor given in components 

by 

1 / Seal ^ 

n — 2\ 2[n — 1) / 

The Weyl tensor of (M, g) is the 4-tensor 

^abcd = -Rabcd + Pac9bd + Pbd9ac — PbcQad — Pad9bc, 

and the Cotton tensor is the 3-tensor 

Cabc = DaPbc ~ DbPac- 

If the metric g is replaced by a conformal metric eg, then the Weyl 
tensor transforms as Wcg = cWg. If n = 3 then ly = 0, but one has 

Ccg = Cg. 

A manifold (M, g) is called conformally flat if some conformal man- 
ifold (M, cg) is flat. Any 2-dimensional manifold is locally conformally 
flat. A 3-dimensional manifold is locally conformally flat if and only 
if its Cotton tensor vanishes identically, and a manifold of dimension 
n > 4 is locally conformally flat if and only if the Weyl tensor vanishes 
(see i). 
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